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O ■ ABSTRACT 

C/j Aims. We seek to constrain the formation of the Galactic bulge by means of analysing the detailed chemical composition of a large sample of 
Cu red clump stars in Baade's window. 

Methods. We measure [Fe/H] in a sample of 219 bulge red clump stars from R=20000 resolution spectra obtained with FLAMES/GIRAFFE at 
the VLT, using an automatic procedure, differentially to the metal-rich local reference star /jLeo. For a subsample of 162 stars, we also derive 
[Mg/H] from spectral synthesis around the Mg I triplet at /I 6319 A. 

Results. The Fe and Mg metallicity distributions are both asymmetric, with median values of -l-0. 16 and +0.21 respectively. The iron distribution 
is clearly bimodal, as revealed both by a deconvolution (from observational errors) and a Gaussian decomposition. The decomposition of 
the observed Fe and Mg metallicity distributions into Gaussian components yields two populations of equal sizes (50% each): a metal-poor 
component centred around [Fe/H] = -0.30 and [Mg/H] = -0.06 with a large dispersion and a narrow metal-rich component centred around 
[Fe/H] = -1-0.32 and [Mg/H] = -1-0.35. The metal poor component shows high [Mg/Fe] ratios (around 0.3) whereas stars in the metal rich 
component are found to have near solar ratios. Babusiaux et al. (2010) also find kinematical differences between the two components: the 
metal poor component shows kinematics compatible with an old spheroid whereas the metal rich component is consistent with a population 
supporting a bar. In view of their chemical and kinematical properties, we suggest different formation scenarii for the two populations: a rapid 
formation timescale as an old spheroid for the metal poor component (old bulge) and for the metal rich component, a formation over a longer 

. ^ time scale driven by the evolution of the bar (pseudo-bulge). 

J^ Conclusions. Guided by these results, we build a simple model combining two components: a simple closed box model to predict the metal 
poor population contribution, whereas the metal rich population is modelled using the observed local thin disc metallicity distribution, shifted 
in metallicity. The. pseudo-bulge is compatible with being formed from the inner thin disc, assuming large (but plausible) values of the gradients 
in the early Galactic disc. 

Key words. Galaxy: bulge - Galaxy: formation - Galaxy: abundances - Stars: abundances - Stars: atmosphere 

1. Introduction tablished. With various stellar population characteristics sim- 
ilar to those of the central old spheroids found at the cen- 

The Milky-Way bulge has been the subject of quite intense ^^^ ^f earUer-type galaxies, and others (mostly geometrical 

debates m the community, as its status is not yet fully es- ^^^j kinematical) that rather point towards a very strong in- 

fluence of a bar responsible for a pseudo-bulge, it sits at 

Send offprint reque.sts to: V. Hill the border between these two types of bulges. The forma- 

* Based on ESO-VLT observations during Paris Observatory tion scenarios for bulges can be classified in three differ- 

FLAMES GTO71.B-0196. ent types: (\) init ial collapse of gas at early times (see e.g. 
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Eggen et al.lll962l) : (ii) merging subclumps, either through an 
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early disc evoluti on (INoguchi 



through mergers (lAguerri et al. 



991 llmmeh et alJ l2004. or 
200 Ic IScannapieco & Tissera 



2003t lNakasato & Nomoto 2003^; (iii) secular evolution of the 



Galac t ic disc (Combes & Sanders 1981; Pfen niger & Norman 
1990t iRaha et alJ Il99lb . The merger scenario itself is simi- 



lar to an early collapse scenario from the point of view of 
the formation characteristic timescale (early and fast). It was 
also recently suggested that the Galactic bulge could be the re- 
sult of both formation processes, with an old sph eroid com- 
plemented by a bar-driven pseudo-bulge (Nakasato & Nomoto 
[20031 iKormendv & KennicuttI 120041: lGerhardll2006l) . The rel- 
ative importance of the two processes (or even populations) 
however remains to be established. 

The prese nce of a bar in the i nner Galaxy has been 
suggested by Ide VaucouleursI (1 19641) from gas kinematics 
and confirmed since then by num erous studies includi ng in- 
frared surface brightness map (e.g. Skrutskie et al. 2006) , star- 
counts, red-cl ump distances (e.g. 'Babusiaux & Gilmore 2005t 
Nishivama et al. 2006; Rattenburv et al. 2007), microlensin: 



and stellar kinematics (e.g. IZhao et al.l 1 19941; [Howard et al 



2008ll2009h . The boxy a spect of the bulge , detected in the in- 



frared light profile (e.g. iDwek et al.lll995h . also argues for a 



pseudo-bulge secularly evolved from the galactic disc. 

On the other hand, photometric studies in selected win- 
dows on the Galactic bulge, in the visible and the near 
infrared, soon made apparent that th e stellar populations 
in t he bulge is old and metal-rich ('Ortolani et al.' ' 19951 
Feltzing & Gilmore 2000; Kuijken & Rich 2002; Zoccali et al. 
2003HClarkson et al.l2008h . Spectroscopicallv. lRichI d 1 988 ) ob- 



tained one of the first metallicity distribution in Baade's win- 
do w based on low-resoluti on spectra of M giant s, followed 
by Ibata & Gilmord (Il995allbl) using K giants and lSadleret al 



(Il996h using red clump stars, all finding a large metallicity 
dispersion. High resolution spectr a of a limited number of 



stars (10 to 20) in Baade's w indow dMcWilham & Richl[l994 ; 
Fulbright et al.l 20061 2007h. and more re cently by our group 
JZoccaU et al.l l2006t iLecureur et al.l l2007b in a larger sample 
of 50 stars in four windows of the Galactic bulge, showed en- 
hanced [a/Fe], compatible with a fast chemical enrichment of 
the Galactic bulge. Both the chemical and age properties of 
stellar popultions in the Galactic bulge thus point towards a 
rapid bulge formation. Early combined metallicity and kine- 
matics measurements also poi nted toward s bulge formation 
through dissipational collapse (iMinnitill 19961) 



More recently still, our group presented in IZoccali et al 



( 120081) the first metallicity distribution entirely based on high- 
resolution spectra of large sample of red giant branch stars in 
three fields of the Galactic bulge (close to the minor axis, at 
b~ -4, -6 and -12°). In this paper we will show how determin- 
ing metallicity distributions from a large (> 200) and almost 
uncontaminated sample of red clump stars in Baade's window, 
observed at high spectral resolution can lead to significant im- 
provements in our understanding of the origin and nature of the 
Galactic bulge. In particular, we show how using two indepen- 
dent elements (iron and magnesium) as metallicity tracers (and 
their ratio [Mg/Fe]) reveals the nature of the stellar population. 
In Sect. 2, we present our target selection and observations with 
FLAMES on VLT, in Sect. 3 we detail the stellar parameter 



and elemental abundance measurement methods, while Sect. 4 
examines the issues of sample representativity and contamina- 
tion. Finally, in Sect. 5, we discuss the resulting [Fe/H] and 
[Mg/H] metallicity distributions and [Mg/Fe] trends, show that 
the sample can be separated into two distinct populations and 
interpret these in the framework of various formation scenarios 
and chemical evolution models. Sect. 6 gives our conclusions 
from this work. 

2. Target selection, observations, and sample 
representativity 

2.1. Target selection 

We selected a sample of stars in the B aade's window field 
BUL-SC45 from the OGLE II survey Judalski et alJ Il997 : 
Paczynski et al.ll999h . among the ~ 1400 stars identified as red 
clump members by the OGLE survey (see Fig. [TJ. The sam- 
ple was first restricted to targets also present in the 2MASS 
dSkrutskie et al.l l2006l) and DENIS surveys dEpchtein et al 



1997h . in order to expand the photometric coverage to the IR, 
which is an asset for photometric temperature determination 
(see Sect. |3]l. From this subsample (~ 800 stars), a further 
restriction was applied in the (Kq, (J - K)()) colour magni- 
tude diagram to lower a possible disc contamination: stars with 
Ko < 12 and (J - K)o < 0.5 were rejected, as they lie in a locus 
where dwarfs are expected to extend. To do so, colours were 

dereddenned using the mean extinction value for the Baade's 
i II 1 

window Ay = 1.5 (Paczynski et al. 1999) and the reddening 
law Ak/Av = 0.089 (Glass. .l999.) . Stars with 16.5 < V < 17.2 
were then selected for observations with FLAMES -GIRAFFE, 
letting the fibre allocation procedure randomly pick the final se- 
lection of 228 stars (see Fig.©- The sample OGLE and 2MASS 
identification and photometry are given in Table [1] available at 
the CDS in its entirety (the first lines of the table are reproduced 
in the printed version). To optimise the exposure time, the to- 
tal sample was divided in two subsamples: (a) 1 14 stars with 
16.5 < V < 16.9 and (b) 114 stars with 16.9 < V < 17.2. In 
each subsample, ~16 fibres were devoted to blank sky regions 
to allow for proper sky- subtraction. In total, 130 fibres were al- 
located in each of the two GIRAFFE fibre configurations (i.e. 
close to the maximum allowed). 

2.2. Observations 

The observations were performed during August 2003 with the 
ESOA^LT/FLAMES facility, as part of the Guaranteed Time 
Observation programs of the Paris Observatory (PI, A.Gomez). 
The spectra were obtained with the GIRAFFE spectrograph 
in Medusa mode using two high-resolution settings: HR13 et 
HR14, yielding resolving powers R = 22500 and 28800, re- 
spectively. Spectral coverage spans from 6120 to 6405 A for 
HR13 and from 6383 to 6626 A for HR14, and the total ex- 
posure times achieved were of 2hl5 and 3hOO in HR13 for 
sample a) and b) respectively, and 4h30 and 6h00 in HR14 for 
sample a) and b) respectively. Using the FLAMES link to the 
UVES spectrograph, a subsample of 12 stars were simultane- 
ously observed at higher resolutions (R ~ 48000). This sub- 
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Table 1. Sample identifications, BVRI OGLE and JHK 2MASS photometry. This table is available electronically in its entirety, 
including also errors on the photometric bands 



ID 


OGLE-ID 


B 


V 


R 


I 


2MASS-ID 


J 


H 


K 






mag 


mag 


mag 


mag 




mag 


mag 


mag 


BWc-1 


BUL-SC45 393125 


20.985 


18.587 


16.838 


15.010 


J18035033+3005324 


13.647 


12.857 


12.677 


BWc-2 


BUL-SC45 545749 


20.434 


18.830 


17.191 


15.392 


J18035671+3005378 


14.115 


13.312 


13.178 


BWc-3 


BUL-SC45 564840 


20.179 


18.689 


16.906 


15.144 


J18035461+3001064 


13.845 


13.087 


12.720 


BWc-4 


BUL-SC45 564857 


19.593 


18.240 


16.760 


15.161 


J18035531+3000576 


13.941 


13.111 


12.886 


BWc-5 


BUL-SC45 575542 


21.116 


18.750 


16.982 


15.175 


J18035592+2955439 


13.909 


13.267 


13.030 


BWc-6 


BUL-SC45 575585 


19.497 


18.238 


16.744 


15.069 


J18035640+2955122 


13.832 


13.097 


12.955 


BWc-7 


BUL-SC45 575585 


19.497 


18.238 


16.744 


15.069 


J18035640+2955122 


13.832 


13.097 


12.955 


BWc-8 


BUL-SC45 78255 


20.702 


18.616 


16.972 


15.163 


J18031236+3003596 


13.724 


12.990 


12.885 


BWc-9 


BUL-SC45 78271 


20.588 


18.527 


16.903 


15.131 


J18031656+3003517 


13.885 


13.153 


12.950 


BWc-10 


BUL-SC45 89589 


19.790 


18.213 


16.695 


15.058 


J18031877+3001101 


13.401 


12.709 


12.505 
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Fig. 1. (I,(V-I)) colour magnitude diagram for the Baade's win- 
dow field BUL-SC45 from the OGLE II survey. The stars iden- 
tified as red clump members by the OGLE survey are indicated 
in blue and the targets observed with FLAMES/GIRAFFE are 
indicated in red. 



sample was previous ly described in IZoccah et al.l ( 120061) and 
Lecureur et al.l (l2007h . 

The data reduction was carried out using the girBLDRS Y_ 
pipeli ne developed at the Geneva Observatory (IBlecha et al. 
2000) which includes cosmic -ray removal, bias subtraction, 
flat-field correction, individual spectral extraction and accurate 
wavelength calibration based on simultaneous calibration ex- 
posures. Co-addition of the individual spectra and sky substrac- 
tion was performed independently from the girBLDRS, using 
several IRAF tasks. For each setup, an "average sky" was made 
by combining the 16 sky fibre spectra and subtracted to the 
spectra of each target. All the spectra for each star were then 
co-added with a /f-sigma clipping to remove the cosmic rays. 



2.3. Contamination and sample representativity 

Since all RGB stars experience the core helium burning phase, 
red clump stars are excellent tracers for the metallicity distri- 
bu tion of a given stellar population. Moreover, as mentioned 
bv lFulbright et al.l ( 120061) . the metallicity-dependent lifetime of 
the horizontal branch (HB) has a negligible effect on the metal- 
licity distribution (MD). These authors estimated that an in- 
crease of 1 dex corresponds to an increase of 10% in the HB 
lifetime, leading to a correction of -0.02/-0.03 on the MD. Red 
clump giants are also good candidates to sample the bulge MD 
because the CMD region where they stand is kno wn to be little 



contaminated b y the other Galactic components dSadler et al 



1996; Fulbrigh t et al.l 120061) . However, since stellar evolution 
models predict that the less metallic stars burn helium in their 
cores on the blue horizontal branch (and not in the clump), it is 
on the other hand expected that there is a minimum metallicity 
that red clump stars can sample. We checked how our sample 
selection centred on the red clump would have missed metal 
poor stars by the in spection of 9 and 12Gyr old isochrones by 



Girardi et al.l (120001) with various metallicities, in J-K, V-I, V-K 



and V-K bands, and compared them with our selection box (al- 
lowing for reddening). We thereby estimated that stars down to 
[Fe/H]> -1.3 or even -1.5 could be included in the sample se- 
lection box. However, because mass-loss is poorly known and 
plays a major role in shaping the horizontal branch, it is dif- 
ficult to quantify, for each metallicity bin, the bias that is in- 
troduced by selecting the red clump. However, the fact that the 
MD decreases smoothly up to our metal-poor limit suggests 
that we might be loosing only a few stars. Furthermore, the 
comparison o f our red clump based MD to that obtained by 
Zoccali et al.l (|2008) from a sample selected higher up on the 



see http ://girbldrs . sourcef orge.net/| 



RGB with no such built-in bias (see Sec. 14.1.41 ). shows no sig- 
nificant difference between the two samples on the metal-poor 

side (that however only reaches to [Fe/H] 1 .2 to - 1 .5 dex). 

We estimated the expected contamination of our bulge 
red clump sample by foreground (and background) Galactic 
disc(s) stars and hal o stars, using the latest population synthe - 
sis Besangon model (Robin et al. 2 0031 : IPicaud & Robin 2004). 
A 25' X 25' field, centred on Baade's Window, was simulated 
without including the reddening. The latter was added indepen- 
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dently of the Besangon model with the following prescription 
based on best reproducing the locus of the CMD in this region 
(position of the foreground disc main-sequence and red clump 
sequences): for stars located beyond 2 kpc from the Sun (in- 
cluding the bulge stars) we applied the same reddening as the 
one used to compute the photometric temperature (Sect. [3] i.e. 
E(B - V) = 0.55), while for stars within 2 kpc from the Sun 
(mostly thin disc), a reddening proportional to the distance (to 
the Sun) was applied (Ay - 0.85/kpc). Further shifts of the 
order of 0.1-0.2 mag had to be applied to both V and V-I in 
order to well reproduce the position of the bulge clump in the 
observed CMD. As shown in Fig.|2] the simulated CMD shows 
a global morphology in agreement with the observed one al- 
though the colour distribution of the simulated CMD is nar- 
rower than the observed, owing to the combined effects of our 
ignoring the differential reddening within the field, and neglect- 
ing the photometric errors in the simulation. 

In the sample selection box (16.5 < V < 17.2 et 1.4 < 
V - I < 2.2), the total number of simulated stars is of ~ 31200 
stars, among which the vast majority are bulge stars (28383). 
The resulting contamination predicted by the Besan9on model 
forour sample is of 9.3%: 5.9% from the thick disc, 3.3% from 
the thin disc and 0.1% from the halo. For the sample of 219 
bulge clump stars, this corresponds to 7+3 stars of the thin disc 
and 13+4 stars of thick disc. It can be appreciated from these 
low numbers that our red clump sample should be very little 
contaminated by any other intervening Galactic populations. 

As illustrated in Fig.|2] the thick disc stars that enter our se- 
lection box are mainly giants, with temperatures and gravities 
similar to those of the red clump stars, and their large majority 
are located beyond 7kpc from the Sun (i.e. in the inner parts of 
the Galaxy). For these thick disc contaminating stars, the anal- 
ysis method developed for the red clump sample and in par- 
ticular the photometric gravity computed assuming that all the 
stars are located at a distance of 8 kpc, can also be applied and 
will provide an unbiased [Fe/H] estimation for these stars. The 
thin disc stars on the other hand, are also giants, but located in 
majority 3-6 kpc from the Sun, and therefore have significantly 
higher gravities than those of the red clump stars. This differ- 
ence in log g reaches 0.4 to 1 dex, which in turn leads to biases 
in the derived [Fe/H] (with our method assuming that all stars 
are located in the bulge, see Sec. [3]l of approximately -0.2 to 
-0.25 dex. Whereas thick disc contaminants from our sample 
could be readily identifiable from their metallicity alone (as- 
suming we know the metallicity distribution of the thick disc 
in the inner regions of the galaxy), thin disc contaminants are 
not expected to be easily identifiable from their position in the 
metallicity distribution. 

The Besangon model makes a number of hypothesis how- 
ever, which may impact the estimates of contamination that we 
make here. One of them is the fact that the thin disc density is 
assumed to decrease sharply in its central 2 kpc (it has a hole in 
the centre), which has been adjusted (together with its scale- 
length) to reproduce star counts to wards the bulge f r om th e 
DENIS experiment, as explained in IPicaud & RobinI (120041) . 
This is the main reason for the little contamination by the thin 
disc predicted here: would the thin disc continue with an ex- 
ponential profile all the way to the galactic centre, the RGB 



and RC stars in the central l-2kpc of the galaxy would con- 
tribute this contamination. On the other hand, observational ev- 
idence of this hole is quite compelling, from star counts, from 
the observation of a such a void in the gas surrounded by a 
molecula r gas ring, and the compelUng evidence for a stel- 
lar bar ((Babusiaux & G ilmord 120051: iNishivama et al.l 12006 : 
Rattenburv et al. 2007)) that would have swept much of the 
matter in these inner regions. More importantly perhaps, con- 
straints on the density profile of the inner parts of the thick disc 
are very poor: the Besangon model assumes a simple exponen- 
tial profile, with no hole in the centre. This results in the major 
part of the predicted contamination to our sample to originate 
in the inner regions of the thick disc: this prediction is therefore 
rather poorly constrained. 

Metallicity distributions assumed in the model also affect 
the position of the stars in each population in the colour mag- 
nitude diagram, and can shift stars in or out of our observed 
selection box. The thin disc is assumed to host a metallicity gra- 
dient of -0.07dex.kpc"' (as constrained from the young stars 
and HII regions), which brings the innermost parts of the disk 
(around 2kpc from the centre, or 6kpc from the sun) to super- 
solar metallicities. If the inner thin disc had no hole, thin disc 
RGB and RC stars in the central l-2kpc of the galaxy would 
have colours and magnitudes (and metallicities) very close to 
our sample's. However, the main-sequence counterpart of this 
younger population have so far elu ded deep colour magnitude 
diagrams of low -reddening regions dClarkson et al.ll2008i) . The 
thick disc contribution is again the most uncertain part of the 
Besangon model prediction, as there are currently no observa- 
tional constraints regarding its radial metallicity gradient (or 
absence of). Would the thick disc in the model have a radial 
gradient, it may bring a few more stars in the selection box. 
However, even with a metallicity gradient, it is unlikely that 
thick disc stars will contaminate the high metallicity end of the 
sample, but would rather be found around metallicities of -0.5 
or below. 

3. Determination of stellar parameters 

3. 1 . Photometric parameters 

Photometric temperatures were determin ed from the V-K, V-H 
and V -J colours using the calibrations of iRamirez & Melendea 
(12005') and their extinction law with E(B - V) = 0.55 to correct 
for reddening. 

The photometric gravities were calculated from the follow- 
ing equation: 

'°^ (i) = '°^ (^) - "" <""■" - "■*■' * "'"^ (^)<" 

with Mbo1,0 = 4.72 mag, Teff.o = 5770 K, loggo = 4.44 dex 
for the Sun and M = O.8M0 for the bulge stars. The bolometric 
magnitude was computed from the V magnitude and the bolo- 
metric correction calibrations of lAlonso et al.l(ll999h assu ming 
that a ll stars were bulge members, at 8 kpc from the Sun (IReidl 
1993h . 

Despite the use of infrared bands, both photometric Tefi and 
log g values are affected by the presence of the differential red- 
dening but show a different sensitivity. In the Baade's Window, 
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Fig. 2. The two upper panels show the (V,(V-I)) colour- 
magnitude diagrams of our field in Baade's window, as ob- 
served (left) and simulated by the Besan§on model (right). The 
red boxes in these panels show the sample selection. On the 
lower part of the figure, the distributions in gravities (left) and 
distances (right) of model stars falling in the selection box are 
presented, separated into their parent populations: bulge (red), 
thin disc (black) and thick disc (blue). The low level of ex- 
pected contamination of our bulge red clump sample is evi- 
dent on these distributions. Note also that, while the thin disc 
contaminants lie in the foreground, thick disc contaminants are 
located in the very inner regions of the Galaxy, at distances 
similar to the bulge red clump population. 



the reddening spread is of the order of 0.15 mag on V-I lead- 
ing to ~200K uncertainties on Tetf and ~0.05 dex on log^. 
Therefore the photometric temperature (mean value from the 
V-K, V-H and V-J colours) was only used as an initial value in 
the stellar parameter determination procedure (see Sect. 13.3b . 
The main source of uncertainty in the log g determination is 
the distance of the star Assuming that all sample stars are at 
the mean bulge distance to the Sun can lead to an error up to 
0.25 dex on \ogg. However, on the whole sample, the previ- 
ous uncertainty value remain lower than those expected from a 
log g deduced from the ionisation equilibrium. This is mainly 
due to uncertainties on the Fe II abundances coming from un- 
certainties in the equivalent widths measurements (continuum 
plac ement, blends) and on the log gf values of the Fe II lines 
(see iLecureur et al.l (l2007h '). At GIRAFFE resolution, the pre- 
vious source of uncertainties become larger and the photomet- 
ric gravity was used as a final value. We also note that the Teff 
value influences the photometric log g: a change of ~ 200 K 
can lead to changes on logg up to 0.15 dex for the cooler stars. 
Then, the log g value were recomputed at each step of the stel- 
lar parameters determination procedure taking into account the 
change on Tgff . 



3.2. Spectroscopic parameters 

3.2.1 . EWs measurement, atmospheric models and 
codes 

The equivalent widths (EWs) for selected Fel and Mgl 
hnes were measured usin g the automatic code DAOSPEQ^ 
dStetson & Pancind 120081) . The synthetic spectra were com- 
puted using t he LTE spectral anal ysis code "turbospectrum" 
(described in I Alvarez & Plea 1 1 998 ) and the abundances from 
EWs were derived using the Spite programs (iSpitdl 19671 and 
subsequent improvement s over the years) with f or both codes 
the new MARCS models JGustafsson et al.ll2008l ). 



3.2.2. The iron linelist 



The iron linelist from ILecureur et al.l (12007ft was used for the 
stellar parameter determination. This linelist was established in 
the following way: a s a starting point, we used the linelist from 



g ' 

Zoccali et a l.^ ('20041 ) with log ;?/ value s extracted from the 
NIST databasqj (R alchenko et alj|2008h . By computing syn- 



thetic spectra in the stellar parameters range of the red clump 
stars, we rejected the iron lines blended up to 10% of the equiv- 
alent width. Finally, for the remaining iron lines, the log gf 
values were adjusted in order that each line gives an abundance 
of 0.30 dex from the EW measured on the observed spectrum 
of //Leo with the following stellar parameters: Teff = 4540 K, 
\ogg - 2.3 dex, and ^ = 1.3 kms '. The spectrum of yuLeo 
was obtained at the Canada-France-Hawaii Telescope with the 
ESPaDOnS spectropolarimeter (resolution R~80000 and S/N 
per pixel ~ 500) an d processed using t he "Libre ESpRIT" data 
reduction package (iDonati et alJll997n . 

Restricted to the wavelength ranges of H13-HH14, the 
linelist contains 48 Fel lines (92 initially). We checked the 
consistency of this reduced linelist on Arcturus and the Sun 
but also on the 12 red clump stars observed with GIRAFFE 
and UVES. The same parameters were found for Arcturus with 
the reduced and total linelists : Teff = 4300 K, logg = 2.5 dex, 
^ = 1.5 kms ' and [Fe/H] = -0.52 dex. For the Sun, we ob- 
tained logg = 4.4 dex, ^ = 1.5 kms"' but Tefr= 5877 K, ie in- 
creased by 100 K to fullfill the excitation equilibrium criterium. 
However, this has no impact on the deduced metallicity. For the 
12 red clump stars, a new determination of the stellar parame- 
ters has been performed with the reduced linelist from the EWs 
measured on the UVES spectra with the same method as de- 
scribed in ILecureur et al.l (120071) . On the mean, the difference 
on the Teff values are null with a dispersion of 100 K, lower than 
the uncertainties. The ^ resulting from the use of the reduced 
linelist are slightly lower (0. 1 km s ') than the ones found with 
the total linelist which translates into [Fe/H] values higher (0. 1 
dex) than the previous ones. However, these systematics are 
smaller than the uncertainties found on the individual values 
and let us conclude that the stellar parameters determination 
from the reduced line is reliable. 



^ DAOSPEC has been written by P.B. Stetson for the Dominion 
Astrophysical Observatory of the Herzberg Institute of Astrophysics, 
National Research Council, Canada. 

^ Available: http://physics.nist.gov/asd3 j 
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Fig. 3. Results of the simulations for the star BWc-2 (Teir = 
4656 K,logg = 2.30 dex,^= 1.5 km s' and[FeI/H] =0.03). 
Left panel: Histogram of ^ values found with the observed EWs 
(shaded histogram) and theoretical EWs (black line). Right 
panel: Histogram of the mean [Fe I/H] abundances. The cor- 
relation between uncertainties in observed EWs and the Fel 
abundances leads to a systematic overestimate of ^ determined 
from the abundances versus observed EWs plot. This overes- 
timate is not significant anymore with the use of theoretical 
EWs. 

3.2.3. The microturbulence velocity determination 



In lLecureur et al.l (120071) . ^ was determined so that lines of dif- 
ferent observed EWs give the same iron abundance. However, 
the observed EW values are sensitive to uncertainties due to 
either the quality of the spectra (S/N, resolving power) or the 
analysis method (continuum placement, FWHM fittings, ... ) 
which translates to correlated uncertainties on the correspond- 
ing [Fei/H]. This correlation between [Fei/H] and observed 
EWs unc ertaintie s leads to a bias towards higher ^ as men- 
tioned bv lMagainl(ll984 . 

In order to estimate this eff'ect on the GIRAFFE sample, 
we performed simulations of the observed EWs and measured 
f for each simulation. This procedure was applied to the star 
BWc-2 (a star representing well our sample: Teu = 4656 K, 
logg = 2.30dex,^= 1.5 kms"' and [Fe I/H] = 0.03). We sim- 
ulated 1000 stars by randomly drawing a set of Fe I lines EWs 
(EWSiand) around the expected EWs (EWsexp) with a standard 
deviation cr - i5ews (given by DAOSPEC). For each 1000 sim- 
ulated stars, Fel abundances ([Fei/H]rand) and corresponding 
uncertainties were computed from the random EWs using the 
stellar atmosphere model adopted for BWc-2. And ^ was modi- 
fied to minimise the slope p of [Fe i/H]i-and versus EWsrand such 
that |p| < (Tp, with cTp the uncertainty on p. The simulation re- 
sults are shown in Fig. [3] The systematic eff'ect on the determi- 
nation of ^ using the observed EWs is clearly illustrated by the 
f distribution histogram. In the mean, ^ has to be increased by 
0.10 + 0. 11 km s"' leading to a corresponding systematic effect 
on the mean [Fe I/H] abundance which decreases by 0.06+0.08 
dex. 

This systematic eff'ect can be suppressed using an abscissa 
linked to the o bserved EW but not aff'ected by random errors. 
MagainI (11984 ) suggested to use the expected EWs (computed 
from the adopted stellar atmophere model and some assumed 
values of ^) and proposed a scheme where the slope of [Fe i/H] 
versus EWSexp is minimised to deduce ^. Such ^ values are then 
free from bias EW abundance error correlations. This scheme 
was applied to the previously simulated sample. As clearly il- 
lustrated by the Fig. |3] the systematic effect on ^ is not signif- 



icant anymore: the mean ^ value is of 1 .50 + 0.06 km s ' and 
the corresponding mean [Fe I/H] value is of 0.02 + 0.03 dex. 

However, this method can introduce another type of bias 
on ^ (and thus on [Fe/H]) if the EWsexp are computed from an 
atmosphere model that is far from reality. In particular, [Fe/H] 
needs to be known to ~ 0.07 dex. For this reason, [Fe/H] was 
estimated as well as possible before ^ was finally constrained. 

3.3. Final stellar parameters determination 

Starting with the photometric temperature, the final Teu was de- 
termined iteratively by requiring the excitation equilibrium on 
Fe I lines (no variation of [Fe i/H] with the excitation potential 
of the line;^'ex- ^ was constrained requiring no trend in the Fe I 
abundance as a function of the expected EW of the lines. At 
each step of the iteration on ^, we checked that the metallicity 
of the atmophere model was the same as the one derived from 
the average Fe I abundance. As explained in Sect. 13.11 the pho- 
tometric gravity was adopted as a final value and recomputed 
at each step of the iteration taking into account the change in 
Teff and [Fe/H]. Stellar parameters and [Fe/H] values were fi- 
nally obtained for 219 red clump stars of the initial sample (228 
stars), and are reported in Table|2l available at the CDS (the first 
lines of the table are reproduced in the printed version). 

This metho d was also applied to the UVES sample of 
Lecureur et al.l (2007), using their full linelist, leading to new 
stellar parameters for each star of the sample. Fig.|4]shows the 
comparison between the new stellar parameters obtained with 
the new ^ determination method and the "old" ones. There is 
no significant change in Teif values, the difference between old 
and new values is of -34 + 1 19 K. As expected with this new 
method, ^ values are on the mean found systematically lower of 
0.12+0. 14 km s"' which translates into [Fe/H] values systemat- 
ically higher of 0.10+0. 09 dex. Compared to the whole sample, 
the red clump stars are more sensitive to the change of ^ deter- 
mination with a mean difference between old and new values 
which reaches 0.25 km s ' . This can be explained by lower S/N 
for these stars compared to the total sample and consequently 
higher uncertainties on the EWs and corresponding Fe I abun- 
dances making these stars more sensitive to the bias affecting 
the ^ determination from the observed EWs. 



3.4. The Stellar parameters uncertainties 

For each star of the red clump sample, an estimate of the uncer- 
tainty on [Fe/H] has been computed from the uncertainties on 
the stellar parameters. The uncertainties on Tgir and ^ have been 
evaluated from the uncertainties on the Fel line EWs given 
by DAOSPEC which translate into uncertainties in Fe I abun- 
dances derived from the EWs. These individual uncertainties 
on the Fel lines were considered in the least squares line fit 
to the points in the (x^x, [Fe I/H]) plane to fix Teff and in the 
(EWs, [Fe I/H]) plane to fix ^ and are contained in the error on 
the line slope. 

Keeping a fixed ^, Teff upper and lower allowed bound- 
aries were estimated by letting the slope p reaches its 1 cr un- 
certainty: p + CTp. The thus determined boundaries are called 
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Table 2. Adopted stellar parameters and measured [Fe/H] and [Mg/H] in the sample stars. 



ID 


Teff 


log 5 


^ 


[Fe/H] 


O-FeClow) 


a-Fe(up) 


[Mg/H] 


crMg(low) 


crMg(up) 




K 


dex 


kms"' 


dex 


dex 


dex 


dex 


dex 


dex 


BWc-1 


4528 


2.09 


1.4 


0.29 


0.14 


0.13 


0.34 


0.08 


0.08 


BWc-2 


4556 


2.25 


1.4 


0.04 


0.10 


0.21 


0.29 


0.10 


0.10 


BWc-3 


4713 


2.23 


1.5 


0.33 


0.17 


0.28 


0.42 


0.15 


0.15 


BWc-4 


4912 


2.28 


1.4 


0.01 


0.10 


0.19 


0.26 


0.10 


0.10 


BWc-5 


4636 


2.21 


1.5 


0.58 


0.12 


0.33 


0.63 


0.14 


0.14 


BWc-6 


4769 


2.19 


1.2 


-0.17 


0.17 


0.30 


-0.00 


0.19 


0.19 


BWc-7 


4769 


2.19 


1.2 


-0.17 


0.17 


0.30 


-0.00 


0.19 


0.19 


BWc-8 


4810 


2.31 


1.4 


0.49 


0.15 


0.17 


0.47 


0.09 


0.09 


BWc-9 


4685 


2.21 


1.5 


0.17 


0.12 


0.24 


0.21 


0.12 


0.12 


BWc-10 


4690 


2.13 


1.5 


-0.13 


0.14 


0.24 


0.20 


0.11 


0.11 
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Fig. 4. Comparison of the stellar parameters ([Fe/H], ^ and 
Teff) obtained with t he new ^ determinatio n method (NEW) 
and those published in Lecureur et al.l (12007 ') (OLD). The black 
dots represent the RGB stars and the open circle the red clump 
stars. The black line shows the mean difference and the dotted 
lines the +lcr dispersion. With the new method, the ^ values 
are systematically lower which implies [Fe/H] systematically 
higher. 



Teff.up and Teffjo and the corresponding uncertainties on Teff 
o"up(Teff) and (crio(Teff) respectively. Note that these temper- 
ature uncertainties need not to be symmetric. In fact, for the 
219 red clump stars, we found: crsup(Teff)= 180 + 90 K and 
o"inf (Teff)= 220 + 90 K and half of the sample shows asymmet- 
ric uncertainties on Teff which can reach 200 K. Lower (crio(^)) 
and upper (crup(^)) uncertainty values on ^ have been estimated 
with the same principle fixing the Tgff value. Contrary to the 
Teff's, these uncertainties are quasi symmetric with |crinf(^)| = 



A modification in Teff (or in ^) value induces a change in 
[Fe/H] which varies from star to star depending on the initial 
stellar parameters but also on the quality of the EW measure- 
ment (S/N ratio, continuum placement, ... ). As illustrated by 
Fig. |5] changing Teff by 200 K implies a change of [Fe/H] 
around 0.10 dex for the coolest and more metallic stars of the 
sample and around 0.20 dex for the hotter and less metallic 
stars. Changes of ^ by 0.2 kms"' lead to similar variations in the 
[Fe/H] value: around O.IO dex for the coolest and more metal- 
lic stars of the sample and around 0.16 dex for the hotter and 
less metallic stars. Overall in our 219 clump stars sample, the 
uncertainties on [Fe/H] associated to Teff uncertainties are of 
0.15 + 0.07 dex on the mean, and those associated to ^ uncer- 
tainties, of 0.16 + 0.06 dex. The individual values are reported 
in the Table[3] available at the CDS in its entirety (the first lines 
of the table are reproduced in the printed version). 

The uncertainty associate d with [Fe/H] were co mputed 
with the formalism described in llVTc William et al.l(ll995i) taking 
into account the errors on the model atmosphere stellar parame- 
ters Teff and ^ and the errors on the Fe I line EW measurements. 
Errors on [Fe/H] coming from the log g value or the metallicity 
adopted to compute the stellar atmosphere model were found to 
be negligible compared to the other sources. Finally, the vari- 
ance on [Fe/H], cr^(Fe), can be written as: 



^'(Fc) = \^\ < + |^ 



dW 



dT 



o-f 



+2 



de_ 
dT 



CT,e, 



(2) 



o-sup(f)l = 0.2 1+0.08 km s- 



with e the Fe abundance, cr}.,, crl and crl the variance in 
the EWs measurement, Teff and ^, respectively and (Tt;^^, the 
covariance between Teff and ^. The covariances between the 
EW measurements and the stellar parameters ^ or Teff are null 

because these variables are independent. The terms (^j cr^ 

and (ll j cri were computed as explained in the previous para- 
graph. The variance in the EW measurements and the associ- 
ated uncertainty on [Fe/H] were estimated from the line to line 
dispersion divided by the number of lines. The covariance be- 
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Fig. 5. Variation in the resulting [Fe/H] with increasing Tetf by 
200 K (lower panels) and ^ by 0.2 kms"' (upper panels) as a 
function of the initial Teff (left panels) and of the initial [Fe/H] 
value (right panels). The hotter and less metallic stars are the 
most sensitive to a modification of Tetf whereas the cooler and 
more metallic stars are the most sensitive to a modification of 



tween Teff and ^ were approximated by measuring the variation 
in the slope fixing Teff (resp. f) due to changes in ^ (resp. Tgff). 

The total upper and lower uncertainty values on [Fe/H] and 
[Mg/H] computed from equation |2] are reported in Table |2] 
while Table [3] details the different contributions of the stellar 
parameter uncertainties (cr^ and ctj) propagated onto the abun- 
dances (cr^(Fe)). On the whole red clump sample (219 stars), 
and as illustrated in Fig.|6l the lower values range from 0.04 to 
0.35 dex with a median value of 0. 17 + 0.05 and show no trend 
with metallicity. The upper values are systematically higher 
with a median value of 0.24 + 0.08 and show an increase with 
the metallicity reaching 0.5 dex for supersolar metallicity stars. 

We also mention that a mean error in [Fe/H] was also es- 
timated by Zoccali et al. (2008) by comparing the measured 
metalhcities in the 5 stars (RGB and red clump stars of the 
Lecureur et alJ (1200 7) 's sample) observed both with GIRAFFE 
and UVES. They found no systematic offset on the [Fe/H] val- 
ues with a scatter of 0.16 dex, a value in agreement with the 
previous estimates from the total GIRAFFE red clump sample. 



3.5. Magnesium abundances 

The Mg abundances were determined from the 6319A triplet 
using spectral synthesis fitting. The synthetic spectra were 
computed using the LTE spectral analysi s code "turbospec- 

and th e mole cu- 
(l2007h. In 



19981) 



trum" (described in lAlvarez & Plez 

lar and atomic linelists described : 

short, we recall that the atomic linelist is the one of the VALD 






Lecureur et al 



database (IKupka et al.lll999l and references therein) with (i) 
solar astrophysical log gf -values for the three lines of the Mg 
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Fig. 6. Left panels: upper and lower errors on [Fe/H] as a func- 
tion of [Fe/H] for the 219 red clump stars. The upper values 
increase with the metallicity with mean values around 0.3 dex 
for the most metallic stars. The lower values show no trend with 
[Fe/H] with a median value of 0.17 dex for the whole sample. 
Right panels: upper and lower errors on [Mg/H] as a function 
of [Mg/H] for the 162 red clump stars selected from the total 
sample. At low [Mg/H], the abundance errors are dominated by 
the line measurement rather than stellar parameter uncertainties 
because 

the Mg triplet lines become very weak. 



triplet and (ii) the A 6318.1 A Cal autoionisation line broad- 
ening adjusted in order to reproduce the Ca I dip in the Sun, 
Arcturus and //Leo. The molecular line lists included in our 
syntheses are the following: C2 ('^C'^C and '^C'^C) Swan sys- 
tem (A-X); CN ('^C'-'N an d i^C'^^N) red system (A-X); TiO y 
and Y systems |Pleg (Il998h . In and near the triplet, the log gf- 
values of some CN lines were modified in order to obtain the 
same Mg abundance from the three lines of the triplet in the 
three reference stars. 



3.5.1 . In the UVES stars: new determination 

Adopting the new stellar parameters and metalhcities (see 
Sect. I3.3l l, the Mg ab undances were determi ned for the UVES 
stars as described in iLecureur et al.l (l2007h with a Ca abun- 
dance computed from the following relation: [Ca/Fe] - 
-0.20[Fe/H] + 0.06. This linear relation was obtained on the 
UVES sample by a least square fit to the [Ca/Fe] values de- 
duced from EWs. 

The mean difference in [Mg/H] between the two determina- 
tions is very small (as expected because [Mg/H] is not sensitive 
to a change of ^): [Mg/H]new - [Mg/H]oid = 0.03 + 0.06 dex 
but becomes more significant on [Mg/Fe] with [Mg/Fe]new - 
[Mg/Fe]oid = -0.06+0.10dex(seeFig|7ll. However, the global 
trend of [Mg/Fe] vs [Fe/H] is almost not affected by this dif- 
ference as illustrated by Fig|7]which shows the new [Mg/Fe] in 
the bulge stars and the [Mg/Fe] found in the thin and thick discs 
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Table 3. Errors in stellar parameters and associated errors on Fe for the sample stars. 



lower 



upper 



ID 


Teff 


logg 


^ 


[Fe/H] 


o-£ 


<T((Fe) 


o-j 


o-T(Fe) 


o-f 


fTfCFe) 


CTj 


o-T(Fe) 




K 


dex 


kms ' 


dex 


kms"' 


dex 


K 


dex 


kms"' 


dex 


K 


dex 


BWc-1 


4528 


2.09 


1.4 


+0.29 


-0.10 


+0.07 


+ 150 


+0.10 


+0.20 


-0.16 


-200 


-0.10 


BWc-2 


4556 


2.25 


1.4 


+0.04 


-0.20 


+0.19 


+ 100 


+0.08 


+0.10 


-0.07 


-150 


-0.07 


BWc-3 


4713 


2.23 


1.5 


+0.33 


-0.20 


+0.16 


+300 


+0.23 


+0.30 


-0.22 


-50 


-0.03 


BWc-4 


4912 


2.28 


1.4 


+0.01 


-0.20 


+0.19 


+ 100 


+0.10 


+0.10 


-0.06 


-100 


-0.09 


BWc-5 


4636 


2.21 


1.5 


+0.58 


-0.30 


+0.30 


+300 


+0.17 


+0.10 


-0.12 


-400 


-0.16 


BWc-6 


4769 


2.19 


1.2 


-0.17 


-0.20 


+0.15 


+300 


+0.26 


+0.20 


-0.13 


-200 


-0.15 


BWc-7 


4769 


2.19 


1.2 


-0.17 


-0.20 


+0.15 


+300 


+0.26 


+0.20 


-0.13 


-200 


-0.15 


BWc-8 


4810 


2.31 


1.4 


+0.49 


-0.10 


+0.10 


+ 150 


+0.13 


+0.10 


-0.07 


-200 


-0.14 


BWc-9 


4685 


2.21 


1.5 


+0.17 


-0.20 


+0.20 


+ 150 


+0.13 


+0.20 


-0.14 


-200 


-0.15 


BWc-10 


4690 


2.13 


1.5 


-0.13 


-0.20 


+0.15 


+200 


+0.18 


+0.30 


-0.21 


-200 


-0.14 



stars bv lBensbv et"aD (l2004Eo05l) and lReddvet"aD(l2006h . In 
the lower metallicities range ([Fe/H]< -0.4), the bulge stars 
are not as separated fr om the thick disc stars than they were in 
JLecureur et al.ll2007l see Fig. 6) with [Mg/Fe] values in bulge 
stars similar to the highest values found in the thick disc stars. 
For -0.4 < [Fe/H] < 0.1, the bulge still shows higher [Mg/Fe] 
values than those found in both discs. At higher metallicities, 
the dispersion has slightly decreased (this comes mainly from 
our use of the linear relation to compute the Ca abundance 
rather than relying on individual -more uncertain- values for 
each star) but remains high, with [Mg/Fe] values from thin disc 
values ([Mg/Fe] ~ .1) to 0.55. The conclusions drawn by 
Lecureur et alJ (120071) from the [Mg/Fe] values are still valid 
with the [Mg/Fe] computed from the new stellar parameters. 
The abundances of O, Al and Na in the UVES stars will also 
be updated using these new stellar parameters and be published 
and discussed in a future paper. 



3.5.2. In the GIRAFFE stars 

From the total red clump sample (219 stars), we selected stars 
according to their uncertainties on the [Fe/H] values in order 
to exclude stars for which Mg would be measured with too 
large uncertainties (mainly stars with low S/N spectra). Only 
stars with crup(Fe) < 0.30 were kept providing a new sub- 
sample of 162 stars. To compute the synthetic spectrum, the 
abundances of C,N, O and Ca are needed for each star of the 
GIRAFFE sample. As a result from the decrease in spectral 
resolution, the C, N and O indicators measurable in the UVES 
spectra are either nondetectable in the GIRAFFE spectra or 
would lead to large uncertainties in the resulting abundances, 
mainly due to uncertainties in the continuum placement (~ 0.5 
dex on [C/H] and ~ 0.3 dex on [O/H]). Therefore, the synthetic 
spectra were computed with the mean values, [C/Fe]=-0.04 and 
[N/Fe]=+0.43, found in the UVES stai's sample. The O abun- 
dance was deduced from the following linear relation: [0/Fe]=- 
0.56 [Fe/H] +0.22 fitted by least squ ares to the [O/Fe] v alues 
found for the UVES stars sample of iZoccali et al.l ( 120061) . The 
Ca abundance was computed as in Sect. 13. 5 .Tl 



Each observed spectrum was first normalised using the con- 
tinuum found by DAOSPEC. The continuum was then adjusted 
by a visual inspection of a 10 A wavelength interval centred 
on the triplet. In order to well identify the possible absorption 
from the CN lines and/or from the Ca I line, two other synthetic 
spectra were overlaid: (i) one spectral synthesis only with the 
molecular lines and (ii) one spectral synthesis only with the 
Cal. This visual analysis also permits to flag or reject some 
stars from the final Mg measurement for one of the following 
reasons: (i) presence of telluric lines affecting the Mg measure- 
ment, (ii) a too low signal to noise ratio in the Mg region which 
can aff'ect the continuum placement and (iii) a strong disagree- 
ment between observed and synthetic spectra. 

At GIRAFFE resolution, blends between the two redder 
lines of the triplet and between the Mgl line at 6318.72 A and 
the Fe I line at 6318.1 A become more important as illustrated 
by Fig|8] This figure also shows that linelist uncertainties do 
not have a significant impact on the Mg abundance determina- 
tion from UVES spectra, but become larger at GIRAFFE res- 
olution and can affect the Mg value. More specifically, in the 
synthetic spectra, the absorption around 6318.2 A is overesti- 
mated and at GIRAFFE resolution, the left wing of the Mg I 
line at 6318.72 A become more contaminated leading to an 
underestimation of the Mg abundance deduced from this line. 
These effects were measured on //Leo by comparing the ob- 
served spectrum convolved to the GIRAFFE resolution with 
the synthetic spectrum computed with the stel lar parameters 
C, N, O and Ca abundances for //Leo from iLecureur et al 



(120071) and different Mg abundances (see FiglS] right panel). 
The abundance deduced from the two redder lines of the triplet 
is [Mg/H]=0.42, whereas the abundance deduced from the line 
at 6318.72 A is 0.2 dex lower. 

Due to the previous considerations, for each star of the 
subsample, two estimations of the Mg abundance were com- 
puted by minimising the x^ values between normalised ob- 
served and synthetic spectra on a wavelength domain restricted 
(i) to the region covered by the two redder lines ([Mg/H]2Lines) 
and (ii) to the region covered by the three lines of the Mg triplet 
([Mg/H]tiipiet)- On the whole sample, the Mg abundance de- 
duced from the triplet is on the mean 0.05 dex lower than the 
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Fig. 8. Left Panel: Observed spectrum (black points) around the Mg I triplet for the B Wc- 1 red clump star obtained with GIRAFFE 
(lower spectrum) and UVES (upper spectrum). The red lines represent the synthetic spectra computed with [Mg/H]=+0.34 and 
[Mg/H]=+0.38 respectively for GIRAFFE and UVES convolution. Right Panel: Comparison of observed (black points) and 
synthetic spectra (red Unes) computed with [Mg/H]=-0.16 for the red clump star BWc-93 and [Mg/H]=+0.42,+0.32,+0.22 for 
/iLeo. The blue dotted lines and the green thick dashed lines represent the absorption of the Ca 1 line and of the CN lines 
respectively. In BWc-93, the contamination by the Cal and CN hnes is negligible and the agreement between the observed and 
synthetic spectra is found for the three lines of the triplet. In/iLeo, the abundance deduced from the Mg I line at 6318.72 A is 0.2 
dex lower than the one deduced from the two redder lines of the triplet. 



one deduced from the two redder lines with a dispersion of 
0.06 around the mean. For the less metallic stars, the differ- 
ence [Mg/H]2Lines " [Mg/H]ti-ipiet Can be negative and the mean 
difference is null. In this metallicity range, the triplet absorp- 
tion weakens and the abundance deduced from the two redder 
lines becomes more sensitive to uncertainties (S/N, continuum 
placement). For stars with [Fe/H] > -0.4 dex, the difference 
starts to be positive and increases with the metallicity of the 
star to reach 0.15 dex for the more metallic stars. Moreover, at 
supersolar metallicities, the [Mg/Fe] values deduced from the 
triplet are more dispersed than those deduced from the two red- 
der lines. In our red clump GIRAFFE sample, to allow for both 
these considerations, we finally adopted the [Mg/H] value de- 
duced from the triplet for the stars with [Fe/H] < -0.4 dex and 
the value deduced from the two redder lines of the triplet for 
the stars with [Fe/H] > -0.4 dex. 

In order to evaluate the Mg determination for the red 
clump stars, we compared the Mg abundances found from 
the GIRAFFE spectra ([Mg/H]GiR) with those found from the 
UVES spectra ([Mg/H]uvEs) for the stars observed with both 
instruments (see Fig. |9l). To increase the statistics of the com- 
parison, we also compute d with the same method the Mg abun- 
dances in 30 RGB stars of lLecureur et al.l (120071) 's sample from 
th e GIRAFFE spectra with the stellar parameters determined 
bv IZoccali et al.l ( 120081) . For the total sample (44 stars), a sys- 
tematic difference between the two determinations is found: 



[Mg/H]GiR - [Mg/H]uvES = -0.06 ± 0.16 dex which does not 
depend on the metallicity. This difference is small, and can be 
mainly explained by the resolution decrease which affects the 
Mg determination via uncertainties in the continuum placement 
and/or uncertainties in the atomic and molecular linelist that af- 
fect the line fitting differently at different resolutions. Both red 
clump and RGB stars show the same behaviour 

However, in terms of comparison with abundances found in 
the discs, this small difference can become significant. Fig. |9] 
shows the [Mg/Fe] ratios found from GIRAFFE and UVES 
spectra of the same stars. The two sets of points show a simi- 
lar global ti-end, although for -0.3 < [Fe/H] < 0, [Mg/Fe]GiR 
ratios are closer to those found in thick discs stars at the same 
metallicity than the [Mg/Fe]uvES were. At higher metallicity 
([Fe/H] > 0.2), the dispersion on [Mg/Fe]GiR values is consid- 
erably lower than that of [Mg/Fe]uvES- We attribute the larger 
dispersion found from the UVES spectra to the result of the 
lower S/N ratios for these spectra. From both sets of mea- 
surements, the bulge stars have [Mg/Fe] ratios on the mean 
higher than those of the thin disc stars, but half of the stars 
fall within the thin disc trend from the GIRAFFE measure- 
ment, which was not the case from the UVES measurement. 
The [Mg/Fe] trend will be further discussed in Sect. 14.11 for 
the totahed clump sample, but w e would like to note a t this 
point that the previous sample of lLecureur etal J ( l2007b plot- 
ted here was a mix of red clump (13) in Baade's Window and 
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Fig. 7. Upper Panel: Difference on [Mg/Fe] found for the 
UVES red clump (open circle) and RGB stars (full circle) with 
the old and new set of parameters as a function of the new 
metallicities. The black line shows the mean value of this differ- 
ence (-0.06) and the two dotted lines shows the +lcr dispersion 
(cr = 0.10). Lower Panel: [Mg/Fe] in the UVES bulge stars 
(red points, computed with the new stellar parameters) com- 
pared with the [Mg/Fe] in the th in (black squares) and thic k 
(blue tria ngles) discs stars from iBensbv et al.l (120041 120051) : 
Reddv et al. (2006). For metallicities from -0.4 to 0.1, the bulge 
stars show higher [Mg/Fe] values than in the two discs. 
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Fig. 9. Lower Panel: Comparison of the [Mg/Fe] found from 
GIRAFFE spectra (green filled circles) with [Mg/Fe] found 



from UVES spectra for 44 bulge stars of Lecureur et al 



(2007)'s sample (open dotted circles). The thin (black filled- 
squares) an d thick (blue open triang le s) discs abu n dance s 



come from Bensbv et al 



open triang le s) discs apu ndances 
(l2004 l2005h : iReddv et al.1 (l2006l) . 



Upper Panel: Difference on [Mg/Fe] found from GIRAFFE and 
UVES spectra for 44 stars (red clump (open circle) and RGB 
stars (full circle)) observed with both instruments. The black 
line shows the mean value of this difference (-0.06) and the 
two dotted lines shows the +lcr around {o- = 0.16). 



RGB stars (44) coming from four different fields located at 
Galactic latitudes of -3°, -4° (Baade's Window), -6° to -12°. 
Among these fields, some are expected to be more contami- 
nated by foreground Galactic discs (thin and thick), and this 
could partly explain the difference between the sample studied 
in|Lecureur et al. (2007) and the pure Baade's Window sample 
studied in the present paper. 

3.5.3. [Mg/H] and [Fe/H] uncertainties for tine red 
clump sample 

The individual uncertainties on [Mg/H] and [Fe/H] were es- 
timated from uncertainties associated with the x^ procedure 
(o"fit), as well as the uncertainties linked to uncertainties on 
stellar parameters (Teff and ^ computed in Sect. I3.4l i. crpix 
was computed using the 6x^ = 1 contour. The uncertainty on 



[Mg/H] due to the uncertainty on Teff and on ^ were computed 
with the following relations: o-([Mg/H])T = (^t^|^) cr(T) and 
o-([Mg/H])f = (^^^^)cr(^) with o-(T) and cr(^) the uncer- 
tainties on Teff and ^ respectively. ( ^ ) and ( ^ | ■ 



were 



estimated from the [Mg/H] dependance to the Teff and ( vari- 
ations found in the UVES stars (see table 10 of iLecureur et al 



20071) . Errors arising from log g uncertainties (an error of 0.3 



dex on log g corresponds to an error of 0.02 on average) were 
found to be negligible compared to the other sources. 
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4. Metallicity distributions 

4.1. Raw metallicity distributions 

4.1 .1 . [Fe/H] metallicity distribution 

The resulting metallicity distribution (MD) for the 219 bulge 
red clump stars is shown in Fig. [10] The mean [Fe/H] value 
is +0.05 + 0.03 dex and the full range of metallicities spans 
-1.13 < [Fe/H] < 0.71 dex. The distribution is clearly asym- 
metric with a median [Fe/H] value of 0.19 dex and shows a 
large proportion of metal rich stars: 65% with supersolar metal- 
licities and 25% with [Fe/H] > 0.38 dex. For [Fe/H] > 0.50 
dex, the number of stars decreases considerably. 

The MD shows two low-level secondary peaks around 
[Fe/H] = -0.55 and around [Fe/H] = -0.20 dex, both with 
very low statistics. While the first remains whatever the size 
and centre of the histogram bins, the second results from a 
small accumulation of stars in a very narrow metallicity range 
around [Fe/H] - -0.20 and therefore becomes more or less 
pronounced according to binning. We checked that this latter 
peak is not related to the two bulge globular clusters NGC6522 
and NGC6528 that are close to the region: the former ha s an in- 
termediate metallicity close to - 1 (ITerndrup et al.ll 19981). while 
the latter, with a mean metallicity of -0.1 (IZoccali et al.ll2004l) 
is too far from the observed window. 
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Fig. 10. Upper left panel: distribution of [Fe/H] for the 219 
bulge clump stars in the Baade's window (thick histogram) and 
distribution of [Fe/H] for the 162 star subsample for which 
[Mg/H] is measured (hatched histogram). Upper right panel: 
distribution of [Mg/H] for the 162 bulge clump stars in the 
Baade's window. Lower left panel: smoothed distribution of 
[Fe/H] derived from the upper panel, with its associated vari- 
ability band. The initial metal rich side of the distribution 
(dashed line) has been stretched ((full line) to make it consis- 
tent with the error law and to achieve the deconvolution (see 
text). Lower right panel: smoothed distribution of [Mg/H] de- 
rived from the upper panel, with its associated variability band 
(see text). 



Table 4. Statistics of the samples: the size n (number of stars), 
the mean a nd median (in dex ) , the s moothing parameter h ac- 
cording to ISheather & JonesI (Il99ll) and the variability band 
width e. 



Subsample 



median 



[Fe/H] 
[Mg/H] 



219 
162 



0.05+0.03 
0.14+0.02 



0.16+0.04 
0.21+0.03 



0.0881 
0.0871 



0.0605 
0.0707 



In order to investigate the significance of low-level sec- 
ondary peaks, we estimate the probability density function 
(PDF) from the raw data, using a kernel estimator This esti- 
mator is: 



^^'~n^h [ h )' 



(3) 



where (X,) are the observed metallicities, n is the sample size 
and h is the window width, also called the smoothing param- 
eter We chose the kernel function /T to be a Gaussian distri- 
bution. The smoothing p arameter /i in Eg. [3] is es timated using 
the scheme described by Sheather & Jones (1991). The estima- 



tor h, as solution of the Eq. 12 in ISheather & Jones (1991) . is 
given in Table |4] for each subsample. 



Th e variability bands, as described bv lBowman & Azzaliru 
(119971 Chap 2.3), are computed to assess the significance 
of the modes in the derived distributions. In the case of a 
Gaussian kernel function, the variability band width around 

V5 is e w 0.2656 (/in)" ' dBowman & AzzaUni|[l997l) . The 



values are given in Table |4] and they are displayed in the lower 
panels in Fig. [10] More d etails about the smoothing method can 
be found in .Rover et al.l (12007L Sect. 4 and Appendix C). 

As far as the secondary modes observed in the [Fe/H] his- 
togram (Fig. [Tol l are concerned, the smoothed distribution to- 
gether with the variability band indicate that the most defi- 
cient one, around [Fe/H] - -0.55 dex, is significant whereas 
the intermediate one, around [Fe/H] = -0.20 dex, is not and is 
rather part of a transition between the two main modes: around 
[Fe/H] = 0.35 dex and [Fe/H] = -0.55 dex. 

4.1.2. [Mg/H] metallicity distribution 

With the analysis method described in Sect. 13.51 Mg abun- 
dances have been determined for a subsample of 162 red clump 
stars built by excluding stars with high uncertainties in the 
[Fe/H] value (only stars with crup(Fe) < 0.30 were kept). As 
illustrated in Fig. [10] the resulting [Fe/H] distribution has the 
same global shape than the one found before from the total 
sample, as was confirmed by a Kolmogorov-Smirnov test (D 
- 0.0793, p-value - 0.5996). The main difference between the 
two distributions is for [Fe/H] > 0.2 where our selection crite- 
ria has excluded a large number of stars, but the [Fe/H] distri- 
bution for the 162 stars still shows the same sharp decrease of 
the number of stars at very high metallicities. 

The resulting [Mg/H] distribution is shown in Fig. [TO] 
Compared to the [Fe/H] distribution, the [Mg/H] distribu- 
tion is much narrower, ranging from [Mg/H] - -0.7 to 
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[Mg/H] - +0.6. While the sample contains 25% of stars with 
[Fe/H] < -0.2, those stars all have [Mg/Fe] > 0, which 
explains the lack of the metal poor tail of the [Mg/H] dis- 
tribution. The mean and median values are higher in [Mg/H] 
than in [Fe/H]: 0.14 and 0.21 respectively (to be compared to 
0.05 and 0.16). The distribution also shows two peaks: one 
around [Mg/H] = +0.3 and a second less significant around 
[Mg/H] = -0.3. For [MgH] > +0.3, the number of stars 
decreases very abruptly. This characteristic has akeady been 
seen in the [Fe/H] distribution but is even more drastic on the 
[Mg/H] distribution. Moreover, we would expect this effect to 
be larger with the addition of the stars excluded by the selec- 
tion. Indeed, the latter are mainly very high metallicity stars 
([Fe/H] > 0.2) that should follow the same [Mg/Fe] trend as 
the other Fe-rich stars (see FigfT4li which have [Mg/Fe] abun- 
dances from -0. 1 to 0. 1 and therefore they should populate the 
region where 0. 1 < [Mg/H] < 0.3. This effect will be discussed 
further in Sect.|5] 

The smoothed distribution of [Mg/H] does not display any 
significant secondary mode, but the spread of the distribution 
is indeed smaller than the one observed for [Fe/H]. 



4.1.3. [Fe/H] distribution robustness 

To further investigate the robustness of our iron MD against 
abundance undertainties, we rejected stars with high uncertain- 
ties in two different ways, and checked the consistency of the 
resulting MD: /) stars with crio(Fe) < 0.22 were kept; //) stars 
with cr„p(Fe) < 0.28 and [Fe/H] < 0.20 or crup(Fe) < 0.36 
and [Fe/H] > 0.20 were kept, taking into account the increase 
of the crio(Fe) uncertainty with metallicity (see Sect. 13.41 and 
Fig. |6]l. The two resulting distributions were found to be fully 
compatible with the total MD. Furthermore, to detect possi- 
ble bias in the MD arising from the analysis method itself, we 
checked the robustness of the MD upon the difference between 
the spectroscopic and photometric Teff . The resulting MD (ob- 
tained for the stars with a temperature difference < 50 K and 
> 150 K) are also compatible with the total MD. 

However, the decrease in the iron MD at the high- 
metallicity end (highlighted by the small number of stars with 
[Fe/H] > 0.5) is too sharp to be compatible with the expected 
measurement errors (see 13.41 and Fig. |6]l and suggests that a 
bias against the highest [Fe/H] could be present. To show this, 
let us extract the stars from the MD with supersolar metallic- 
ities, as representative of the metal-rich peak of the distribu- 
tion. First, using [Fe/H] > as a lower cut-off, the cumula- 
tive distribution of supersolar metallicities is plotted in Fig.fTTI 
Under the hypothesis that, even if the highest metallicities are 
biased, the median of the distribution of the metal-rich peak 
is conserved (which is the case if errors remain almost sy- 
metrical), this median (i.e. frequency = 0.5), corresponds to 
[Fe/H] - 0.35dex. This median is relatively robust to different 
cut-off values (-0.1 orO.l dex) for defining the high-metallicity 
peak, for which we obtain respectively [Fe/H] = 0.32 and 
0.36 dex. The frequencies 0.16 and 0.84 respectively mark the 
1-cr positions from the median on the lower and upper sides 
of the distribution. These values, displayed in Fig. (TT] show a 



clear asymetry: crieft = 0.21 dex and (Tnght - 0.16 dex. While 
(Tieft is slightly larger than the median expected measurement 
error criow,Fe/H]>o of 0.17, allowing for a small astrophysical 
scatter, the crjight is too small compared to the median expected 
error o-up,[Fe/H]>o - 0.24 (Fig|6|l. The metal rich end of the 
[Fe/H] MD is thus too steep to be consistent with the expected 
dispersion of the measurements due to their observational un- 
certainties only. One possibility is that remaining degeneracies 
between ^ Tefi and [Fe/H] introduce a bias in the metallicity 
estimates of the most metal-rich stars ([Fe/H] > 0.35) at the 
resolution a S/N of the present sample, preventing us to truly 
establish the underlying shape of the iron MD tail at high metal- 
licities. 



This bias is not detected in [Mg/H] though, for which the 
observed [Mg/H] distribution is compatible with the expected 
uncertainties, even if this distribution also displays a sharp cut- 
off at high metallicities. This is because the Mg lines used in 
the analysis are weak and do not suffer in the same way of ^ 
Teff degeneracies in the stellar parameter determination. As a 
result, the uncertainties on [Mg/H] are both smaller than those 
on [Fe/H] at high metallicities, and less assymetric. Another 
way to say this is that, since the error on [Mg/H] is dominated 
by actual line measurement errors (synthetic spectra fits in this 
case), the possible (small) bias in [Mg/H] at the highest metal- 
licities is hidden. 



To understand wether such a bias could indeed be present, 
we computed synthetic spectra of metal-rich (0. to +0.75dex) 
red clump stars, sliced and convolved them to the same wave- 
length domains, resolution and sampling as our observed 
GIRAFFE spectra. We added photon noise to the spectra to 
reproduce a set of spectra with S/N=100, 50, 40 and 30(100 re- 
alisations of the noise were generated for each synthetic spec- 
trum and noise level). We then retrieved the stellar parameters 
and metallicity using the same method as for the bulge stars 
(equivalent width measurements and iteration on the stellar pa- 
rameters). The results of these extensive simulations show that 
for supersolar metallicities, there is indeed a slight bias in the 
retrieved metallicities, although it is in the contrary direction 
than what would be suggested by the error analysis above: in 
the mean, very metal-rich stars are found slightly too metal-rich 
by our method. In the mean, a [Fe/H] =+0.5 star will be found 
+0.14 to +0.17 dex more metal-rich, whereas a [Fe/H]=+0.25 
star will be found +0.10 to 0.14 dex more metal -rich (the bias 
steadily growing with decreasing S/N from 100 to 30). At and 
below solar metallicities, the bias vanishes. Finally, we also 
confirmed this bias by performing the same test on the observed 
//Leo spectrum, degraded in resolution, wavelength coverage 
and photon noise. In this case, the bias is +0.13 at S/N=50 and 
+0.18atS/N=30. 



Based on these findings, we warn the reader that there may 
be a bias in the highest metallicity stars of the sample, although 
this bias is not understood well enough to be corrected. 
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Fig. 11. The cumulative iron metallicity distribution for 
[Fe/H] > (solid line). The median and the left and right stan- 
dard deviations are derived from the frequency values, assum- 
ing Gaussian errors. 
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Fig. 12. Iron metallicity distribution (MD) of RGB stars ob- 
tained with the automatic procedure (red shade d histogram) 
compared with the MD of Zocc ali et al.l (120081) (Left panel, 
black line) and with the MD of the red clump stars {Right panel, 
black line). The new RGB MD is in total agreement with the 
red clump one {Right panel) but is shifted toward h igh metallic- 
ities compared to the MD of IZoccah et al.l (120081) {Left panel). 
This small difference is due to the difference in analysis of the 
two samples as explained in the text. 



with respect to the other (the RGB MD was on average 0.07 
more metal poor than the red clump). 

The difference between the two MD could arise from small 
differences in the analysis itself. Indeed, even if using the same 
criteria to determine the stellar parameters, the analysis of the 
RGB samples is different on two points: i) the photometric tem- 
perature is estimated from the V-1 index or from the TiO index 
(whereas it is computed from V-1, V-J, V-K for the red clump 
sample) and ii) the stellar parameters are determined in a man- 
ual way (whereas the global procedure is automatic for the red 
clump sample). To investigate the two previous differences, we 
re-analysed the RGB sample using the same automatic proce- 
dure with two valu es for the photometr ic temperature: i) the 
same as the one of lZoccaliet al.l (l2008h and ii) a value com- 
puted from the indices V-J, V-H, V-K (as for the red clump 
sample). From the 204 RGB stars sample, we excluded 6 stars 
as members of the globular cluster NGC6522 and 4 stars as 
suspicious binaries and for the test ii), we also excluded stars 
for which the 2MASS survey only gives lower limits in one or 
more IR bands. Finally, [Fe/H] values were obtained for 194 
stars in the case i) and 166 stars in the case ii). Before mak- 
ing the comparison be tween old and new [Fe/H] values, we 
checked that the MD of lZoccali et all (l2008l) reduced to the 166 
stars of ii) was totally compatible with the total one (204 stars). 
As shown by Fig. [12] the MD obtained for the 166 stars from 
the automatic proce dure has a global shape close to the one of 



Zoccah et al.l (120081) (see Fig.[T2ll with many more stars at very 



high metallicities ([Fe/H] > 0.4) leading to mean and median 
values (0.02 and 0.11 respectively) slightly higher than values 
previously found (-0.01 and 0.05 respectively). However, in 
the cases i) et ii), from the results of a KS test, the MD are not 
compatible with the one of lZoccali et alJ (l2008r) : p-value=0.03 
for i) and p-value=0.05 for ii), but they are compatible with the 
red clump one: p-value=0.65 for i) and p-value=0.58 for ii) as 
illustrated by Fig.[T2]in the ca se ii). These results show that the 
difference between the MD of lZoccali et alj (12008 ) and the red 
clump sample comes from the automatic procedure itself rather 
than the initial photometric temperature adopted. 

Note, however, that the difference between the RGB MD 



presented here and that of IZoccali et al.l (l2008h remains very 
small, and that the conclusions of their paper remain fully valid. 



4.1 .4. Comparison with tine [Fe/H] distribution of RGB 
stars 



In IZoccah et alj ( 120081) . another bulge iron MD was obtained 
in Baade's Window from a sample of 204 bulge red giants 
stars. This MD was found in agreement with the red clump 
stars one and the small differences between the two were found 
to be small enough to combine them to establish the Baade's 
Win dow MD for the com parison with the other fields (see Sect. 
6 of Zo ccali et al. (2008)). However, for other statistical works 
(see .Babusiaux et al...201Q) . the RGB and red clump MD had 
to be totally compatible to be combined but a Kolmogorov- 
Smirnov (KS) test did not confirm that the 2 samples could be 
drawn from a parent population ( p-value - 0.0005). At vari- 
ance, the two samples were found to be slightly offsetted one 



4.2. Unravelling two populations 

4.2.1 . Tine deconvolved metallicity distribution 



From the PDFs estimated in Sect. l4.1.Tl and plotted in Fig.[TOl 
one can rectify the error law effect. However, we have argued 
in Sect. 14. 1.3] that the metal rich end of the iron metallicity dis- 
tribution may be biased at the highest metallicities. In order to 
avoid any spurious result in deconvolving the observed PDF, 
we have chosed to strechthe metal-rich side of the MDF, so 
that cTjigijt finally matches crieft. The stretched distribution is su- 
perimposed in Fig. [To] This stretching will not erase any bias 
in the highest metallicities if it is indeed present, but helps the 
deconvolution algorithm to perform nominally. 

We used the Lucy-Richar dson deconvolution technique 
dLucvi 1 19741 : 1 Richardson! 1 19721) to recover the rectified MD. 
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More detail s abou t the deconvolution method can be found in 
Rover et all (l2007l Sect. 4 and Appendix C). 

In this work, the error law is assumed Gaussian. We chose 
to model the error law as a constant, using the median of the 
(o"iow(Fe) + crup(Fe))/2 and o"iow(Mg): 0.21 dex and 0.14dex 
respectively. These values are used to estimate the standard 
deviation in the Gaussian error distribution, during the Lucy- 
Richardson deconvolution. The stopping criterion described by 
Lucvl (119941) was used and the resulting number of iterations 
was 1 1 and 5 for the iron and magnesium distributions respec- 
tively. The resulting PDF are plotted in Fig. [13] The decon- 
volved [Fe/H] distribution shows a clear bimodality with a low 
metallicity peak around -0.45 dex and a high metallicity peak 
around +0.3 dex. As far as the [Mg/H] is concerned, the decon- 
volution does not show clear features in the resulting distribu- 
tion. This is due to the fact that [Mg/H] abundances are less 
spread and that their associated errors are smaller 
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Fig. 13. Upper left panel: Smoothed distribution of [Fe/H], as 
in Fig. [TOl with the resulting deconvolved distribution using a 
constant error (solid line). Upper right panel: Smoothed distri- 
bution of [Mg/H], as in Fig. [TO] with the resulting deconvolved 
distribution using a constant error (solid line). Lower left panel: 
Resulting Gaussian populations in the [Fe/H] distribution. See 
Table|5]for the parameters of the Gaussians. The error bars rep- 
resent the uncertainty in the amplitude of the Gaussian popu- 
lations. Lower right panel: Resulting Gaussian populations in 
the [Mg/H] distribution. 



4.2.2. Gaussian mixture 

The metallicity distributions described in Section 5.1 sug- 
gest that they may be made up of contributions from dif- 
ferent stellar populations. In a first attempt to identify 
possible sub-structures, we have decomposed the observed 
[Fe/H] and [Mg/H] distributions into a finite number of 
Gaussian components. Th e SEMMUL algorithm, developed by 
Celeux & Dieboltl (Il986l) . was appHed. SEMMUL (Stochastic, 



Estimation, Maximisation, MULtidimensional) is an iterative 
method for numerically approximating maximum likelihood 
estimates for incomplete data with a stochastic step to accel- 
erate the convergence. It allows separation of Gaussian com- 
ponents without a priori knowledge of the number of compo- 
nents (only a maximum number is necessary). In addition, it 
does require defining a set of i nitial conditions. The algorithm 
has been modified bv lArenoul (1 1 993h in order to take observa- 
tional errors into account, thus the modified version allows to 
estimate the cosmic dispersion of the population. 

SEMMUL has been applied to the 1 -dimensional [Fe/H] 
and [Mg/H] distributions separately, and to the 2-dimensional 
joint distribution ([Fe/H], [Mg/H]). Thousands of runs were car- 
ried out in each case, assuming that the sample is a mixture of 2 
or 3 discrete components. In all cases, the algorithm converges 
to a two-component solution with stable results, dividing the 
sample in roughly two equal-sized populations. 

The resulting two components (that we call A and B) are 
summarised in Table |5] For each component, the mean chem- 
ical composition, the dispersion and the fraction (in percent) 
of stars in each group, with the corresponding standard er- 
rors, are given. The A and B population characteristics are very 
consistent among the three separations (one-dimensional along 
[Fe/H] or [Mg/H], or 2-dimensional). 

In Sect. 14.1.31 we have argued that the high-metallicity tail 
of the iron distribution may be biased. In this context, we may 
wonder whether the mean values and dispersions obtained in 
the Gaussian decompositions of population B are realistic. The 
comparison between the estimated mean values ([Fe/H] - 0.32 
and 0.29, see Table |5]l and the median values obtained from 
the cumulative metallicity distribution (for [Fe/H] > 0) in 
Sect l4X3l ([Fe/Hl = 0.32 to 0.36) shows that they are con- 
sistent. On the other hand, the obtained cosmic dispersions 
of population B (0.11-0.12 dex) could arguably be underesti- 
mated. However, according to the same cumulative distribution 
presented in Sect. 14.1.31 the total metallicity dispersion is about 
0.2 dex in this population, and the expected measurement errors 
are of the same order, which implies that the cosmic dispersion 
of the population is small, as obtained by our algorithm. 

4.2.3. Ciiaracterising tine two populations 

Two populations are clearly identified by the Gaussian separa- 
tion exercise. The metal-poor component (A) is centred around 
[Fe/H] = -0.30 and [Mg/H] = -0.06, with a large disper- 
sion (~ 0.3 dex), while the metal-rich component (B) has sim- 
ilar mean Fe and Mg content ([Fe/H] = 0.32 and [Mg/H] - 
-1-0.35) and a very small dispersion (0. 11 in [Fe/H] and 0.07 in 
[Mg/H]). 

The dispersion of population B is very small, comparable 
in fact to the metallicity dispersion of Galactic disc in the solar 
neighbourhood (or even smaller). The more metal-poor popu- 
lation A is on the contrary quite extended, both in [Fe/H] and in 
[Mg/H]. This dispersion could in part be due to contamination 
(see Sect. 12.3b by the thick disc in the inner Galactic regions, 
that could account for ~ 10 - 15% of population A (5.9% of 
the total sample). In fact, we will see in the following section 
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that the chemical signature of this population (A) is not distinct 
from that of the thick disc, and we can therefore not exclude 
that part of it is made up of thick disc stars, indistinguishable 
from the bona-fide bulge population by any of our observables. 

It is further interesting to note that the two components thus 
separated are more widely separated in [Fe/H] than in [Mg/H], 
as reflected by the different mean corresponding [Mg/Fe] ratios 
(see Sect. 14. 31 and Fig.O, hinting at two physically separated 
populations. While the metal-poor component A has a clearly 
defined magnesium overabundance of [Mg/Fe] ~ 0.24, the 
metal-rich component B has a significantly lower [Mg/Fe] ~ 
0.05, not significantly different than the solar neighbourhood 
disc magnesium to iron ratio. 

We tentatively explain the existence of t hese bulge pop 



Illations as the result of different origins. In Bab usiaux et al. 



(1201 0). we combine the abundances and the kinematics (radial 
velocities plus proper motions from OGLE-II) and show a sig- 
nificant difference in the kinematics of the metal rich and the 
metal poor components. The velocity ellipsoid of the metal rich 
component shows a vertex deviation consistent with what is ex- 
pected from a population kinematically supporting a bar The 
metal poor component on the other hand shows no vertex devi- 
ation, consistent with an old spheroidal population. We there- 
fore relate the richest population to a bar-driven pseudo-bulge 
and the more metal poor one to an old spheroid with a rapid 
time-scale forma tion. Similar results ha ve been suggested by 



Soto et al.l (120071) and IZhao et al.l (119941) . based on kinematics 



and low resolution metallicities. 

Let us note finally that this mixture of populations in 
Baade's Window is fully compatible with the variation of 
the metalli city distribution as a function of Galactic latitude 
noticed bv IZoccaU et al.1 (l2008l) (from b=-4 to -12°), that 
we would therefore attribute to the gradual disappearance of 
the pseudo-bulge component B as one moves away from the 
Galactic plane, while the old spheroid (A) would dominate 
out to higher l atitudes. The abse nce of gradients in the inner 
bulge found bv lRich etal.l (120071) is also consistent, as this re- 
gion would be dominated by component B only. We also note 
that population B also has to be predominantly old (perhaps 
not surprisingly if dis c formation occurs inside out), given that 



Clarkson et alj (120081) find a pure old population in a field at 



b - -2.65° (i.e. significantly closer to the plane than Baade's 
Window). 

4.3. The [Mg/Fe] trend in the bulge 



Our Mg results in terms of [Mg/Fe] against [Fe/H] and [Mg/H] 
are shown in Fig. [14] (red points). The [Mg/Fe] ratios in the 
bulge decreases with both [Fe/H] or [Mg/H] from ~ H-0.5 to 
near solar values. Comparing our resu lts with the rece i it opti - 
cal and IR stud i es of b ulge giants frorn Rich & Origlial (l2005h . 
iFulbrigh t et~ail (l2007b and Rich et al.1 (l2007h on the common 
metallicity ranges, we note that the general trend is similar, al- 
though our much larger statistics allow us to define the trend 
with more confidence. 

To understand the relation that the bulge bears to the 
Galactic disc(s), it is interesting to compare the [Mg/Fe] ra- 



Table 5. Separation of the sample in two Gaussian components 
(A and B), with [Fe/H], [Mg/H] or both simultaneously as dis- 
criminant variables. (crpe/Hj and cr[Mg/H] are the cosmic disper- 
sion of the population) 
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tios in the bulge sample to the abundance trends found in 
the Galactic thin and thick discs. In Fig. [141 we plotted the 
thin (black symbols) and th ick (blue symbols) di s c stars abun 



dances from the studies of Bensbv et al.l (l2005h : iReddv et al 



(l2006l) and lFuhrmannI (l2008h . These studies are based on very 



careful and thoughrough analysis of large samples of dwarf 
stars in the solar nei ghbourh ood. As recently p ointed out by 
Melendez et"aLl JJOOSh : lAlves-Brito et alj (I2OI0I) . it would be 
better, for the sake of minimising possible systematics in the 
abundance analysis, to compare bulge giants with solar neigh- 
bourhood thin and thick disc giants observed and analysed in 
the same conditions. However, from a statistical point of view, 
the available samples of loc al giants available for su ch com- 



parisons is still quite small ( Alves-Brito et al.ll2010l ~20/30 



thick/thin disc stars respectively, ). We therefore decided to 
stick to a comparison to local disc dwarfs, once a number of 
verifications were performed: (i) we analysed the [Mg/Fe] of 
two reference stars, jt/Leo and Arcturus ([MgfFe]-+Q.l2 and 
-f-0.38), that respectively fall on top of their parent popula 



tionas derived from local dwarf samples; (ii) IFulbright et al 



(■2007) analysed a sample of local K giants, and found no differ- 
ence between the [Mg/Fe] trend deriv ed from these giants and 
those defined from local dwarfs; (Hi) iMishenina et al.l (120061) 
also analysed a large sample of local red clump giants, and 
reach a similar conclusion for [Mg/Fe]. From these checks, 
we conclude that it is fair to compare the [Mg/Fe] abundances 
for the bulge red clump giants presented here to those of lo- 
cal disc(s) dwarfs. In Fig. [141 the two lower panels show that 
indeed our conclusions would not be altered in any \yay if com- 
paring our sample to local giants ( Alves-Brito et alj|2010l) . We 
however caution that in the highest metalUcity regime, system- 
atic differences of the order of 0. 10-0. 15 on the [Mg/Fe] trends 
are found between diff eren t studies of dwarfs stars (for exam- 
ple l^ddxeFaD llooi) vs. lFuhrmannl(l2008l)). be tween dwarf 



and giant stars sample s (for example iBensbv et al.l (l2005h vs. 



Mishenina et al.l (l2006|)^ and between differ en t disc giants stars 



sampl e (for example IFulbright et alJ (l2007h vs IMishenina et al 



(l2006h ). This is clearly a domain where abundances have to be 
taken with caution in general, and even more so when they are 
derived, as here, from moderate S/N and resolution spectra. 
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Fig. 14. [Mg/Fe] as a function of [Fe/H] (right panel) and [Mg/H] (left panel) for our bulge red clump sample (red points) 
compared with thick (bl ue) and thin (black) discs stars. In the two upper pa nels, the d isc reference sam ples are from the main 
sequence stars studies of Bensbv et al.l (120051) (trianglesl lReddv et alJ (120061) (squares ) iFuhrmannI (120081) (crosses), whereas the 
two lower panels show the solar neighborhood disc giants by iAlves-Brito et alj (120101) . as well as their bulge giants (large filled 
cyan circles). In the low metallicity range ([Fe/H] < -0.3 and [Mg/H] < 0), the thick disc and the bulge show similar [Mg/Fe] 
ratios, higher than in thin disc stars. 



When compared to the local Galactic discs abundances 
these various bulge samples show some differences depending 
on the metallicity range: 

In the range {FejH} < -0.3, the number of stars in the 
present sample has increased dramatically compared to our pre- 
vious UVES sample (Fig. |7]i and they clearly show [Mg/Fe] 
ratios similar to those of the thick disc stars for all [Fe/H]. 
In the same meta l hcity range, our results differ from those 
of iFulbright et al.l (l2007b who found higher [Mg/Fe] values 
than ours and a [Mg/Fe] trend higher than that of their sam- 
ple of disc giants. In fact, i f we use Arctu r us as a reference 
to insure that the results of IFulbright et al.l (l2007b are on the 
same scale as the present work, an offset of ~ 0.10 is ex- 
pected ([Mg/Fe]Fuibright - [Mg/Fe]pi-esent = +0.1 1), bringing the 
[Mg/Fe] of the two studies into agreement. 

In the range -0.3 < [Fe/H] < +0.1, the bulge [Mg/Fe] 
trend from our sample is clearly distinct and higher than that 
of the thin disc, whatever comparison sample is used. In the 
same metallicity range, our bulge stars show [Mg/Fe] values 
on the mean higher than those of the thick disc stars. This 
confirms, with better statistics our p revious results from the 
UVES sample (ILecureur et al.ll2007h as well as the results of 



Fulbright et al.l (120071) . Given the residual systematic effect in 
the [Mg/Fe] determination from GIRAFFE and UVES spectra 
for the stars in common (see Fig. |8]l, the difference between 



thick disc and bulge would be even more pronounced if we 
adopted the UVES abundances as the reference. However, this 
difference between bulge and thick disc has to be viewed with 
caution because of the small number of thick disc stars in this 
metallicity range, and the stil l debated thick-disc nature of this 
high-metallicity tail. Indeed, iBensbv et al.l (120051 l2007h argue 
that the thick disc extends at least to solar metallicities and 
shows a clear decrease of [c/Fe] with increasing iron content 
for stars with [Fe/H] > -0.3 denoting an extended star forma- 
tion period (SNIa enrichment). On the contrary, several authors 
(Reddv et a l. 2006; Fuhrmann 2008) argue that the thick disc 
does not extend significantly in this metallicity regime, and that 
there is no evidence of decreasing [a/Fe] in the thick disc at all. 

In the range [Fe/H] > +0.1, the bulge shows [Mg/Fe] ra- 
tios which are similar to those of the local thin disc, solar on 
the mean, and with a decreasing [Mg/Fe] trend with increas- 
ing metallicity. This result confi rms, with better statistics, the 
results found for the oxygen by iMelendez et al.l (120081) in the 
same metallicity regime. On the contrary, our bulge [Mg/Fe] 



ratios are significantly lower than those of IFulbright et al 
(I2OO7I) . even allowing to shift the latter by ~ -0.10 to insure 
that both studies lie on the same scale. We further note that 



FulbrJght et al.1 (l2007h derived [Mg/Fe] ^ +0.02 for ^uLeo, 0.2 
lower than what we found for this star dLecureur et al.ll2007l) . 
At such high metallicities, both Fe and Mg measurements re- 
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main quite uncertain which could ex plain part of the differ - 
ence between our study and the one of iFulbright et alJ OOOTh . 
Moreover, the comparison of the bulge and disc trends should 
be taken with caution in that metallicity regime, since, as 
pointed out above, systematic differences of the order of 0.10- 
0. 15 are found between the [Mg/Fe] trends at supersolar metal- 
licities between different studies. More specific ally, taking as a 
compa rison the local disk sample of giants by IFulbright et al.l 
(l2007h . we would find that the [Mg/Fe] of the bulge lies above 
that of the local disk, at variance with the conclusion that is 
drawn from the comparison to loc al dwarfs. However, w e do 
not favour this interpretation since 'F ulbright et al.l (120071) also 
finds yuLeo +0. 1 1 dex lower in [Mg/Fe] than we do. 

Very recently, several papers have measured detailed abun- 
dances of a few bulge dwarfs from high resolution and high 
signal to noise rati o spectra obtained du r ing a gr a vitational mi- 
crolensing event (Bensbv et al. 2010b, l2009bllal : ICohen et al. 



20081 I2OO9I: Ijohnson et aDr2007t i2008l) . There are in total by 
now 14 microl ensed dwarfs and su bgiants analysed in the 



Galactic bulge jBensbv et al.ll2010bh . and although the early 



(2007 and 2008) microlens events yielded a sample of stars 
almost exclusively me tal-rich that was irreconcilable with the 
MDF of bulge giants dCohen et al.ll2008l l2009t lEpstein et al" 
I2OIQ1) . the 2009 and 2010 events have uncovered a signifi- 
cant metal -poor population, tami ng down this conclusion to the 
point that iBensbv et al.l (1201 Obi) claims that the two distribu- 
tions could statistically arise from the same parent population, 
albeit a rather small probability. The current MDF from these 
microlensed un-evolved stars is now highly bimodal, with a 
peak at [Fe/H] ~ H-0.3 and another one at [Fe/H] ~ -0.5 
dex. This is very similar to the peaks uncovered by our pop- 
ulation separation in the red clump sample. Interestingly, the 
metal-rich (similar to our population A) and metal-poor (sim- 
ilar to our population B) microlensed stars also have differ- 
ent ages, the latter being a clearly old population, while the 
former spans a range of ages (2-13 Gyrs). This is again com- 
patible with our suggestion that population B is related to the 
old bulge while population A is kinematically related to a bar 
and would be made of a pre-existing disc. Even more inter- 
estingly perhaps, we are now in a position to compare the 
[Mg/Fe] trends in the bulge from dwarf and giants stars: as is 
our bulge giants, [Mg/Fe] in dwarfs overlays the local thin disc 
at super-solar metallicities, and the local thick disc at metal- 
licities {FejH} <- -0.3. The current sample of microlensed 
giants only contains two sta rs in the intermediate metallicity 
range dBensbv et al.l uOlObI) . where we find the bulge to be 
richer in [Mg/Fe] than both the local thin and thick discs. These 
two stars seem to trace the upper envelope of the local thick 
disc. To draw any further conclusion on this metallicity range, 
more microlensed dwarfs will need to be analysed. 

To follow up on our tentative sample separation in two dis- 
tinct populations we colour-coded in Fig. [15] the stars belong- 
ing to the populations A and B as separated in Table |5] along 
[Mg/H], using the probability with which a star was assigned to 
its population. All stars with very high probabilities (> 85%) to 
belong to population A are coded in light grey, while those with 
high probability to be in B are coded in black. Stars that are in 
dark grey are those with low probabilities to belong to one or 



the other class. Population A shows high [Mg/Fe] ratios, with a 
mean of [Mg/Fe] - 0.25, while population B has lower Mg en- 
hancements, close to solar, as discussed above. The intermedi- 
ate RC population (those stars that could not be unambiguously 
classified) has both high and low [Mg/Fe] ratios. In fact, the 
[Mg/Fe] histogram of this population shows two well defined 
peaks, one at -1-0.25 and the other at +0.05, coinciding precisely 
with the means for population A and B respectively. It therefore 
seems that this intermediate population is a mix of stars with 
the same chemical properties of A or B population, rather than 
a smooth transition between the two. This is also visible in the 
run of [Mg/Fe] with metallicity (whether [Fe/H] or [Mg/H], 
Fig. [TSl l at intermediate metallicities, where the stars seem to 
cluster around two discrete [Mg/Fe] values extending the trends 
of population A and B, rather than follow a smooth decreasing 
trend. As a result, we argue that population A, that we associate 
with the old Mg-enhanced bulge, extends up to metallicities 
of at least [Fe/H]=-i-0.1 and [Mg/H]=-i-0.35, and possibly up 
to [Fe/H]=-i-0.25 and [Mg/H]=-i-0.5. Conversely, population B 
probably extends down to [Fe/H] ~ [Mg/H] = 0.0, with mod- 
est Mg enhancements. This will have an impact on the chemi- 
cal evolution models that may be appropriate to represent these 
two populations. 

5. Chemical evolution and bulge formation 

In this section, with the aim of understanding the formation his- 
tory and chemical evolution of the galatic bulge, we compared 
our bulge MD with the results of two r ecent bulge formation 
models: the chemical evolution mo del of Ballero et al.l (l2007bl) 
and the chemodynamical model of llmmeli et al.l (120041) . 



5.1 



Comparing to a chemical evolution model of the 
bulge 



Fi gure fT6| upper panel s hows the MD predicted by the model 
of iBallero et al.l ( l2007bl) together with our observed iron MD. 
For the sake of comparison, the model MD has been convolved 
with the mean [Fe/H] uncertainties of the observed sample 
(Fig. |6]l. By constraining their model p arameters frq r n bulge 



observed abundances ratios and Fe MD, IBallero et alj (l2007bl) 
concluded that the bulge formed in a very short time, with a 
high star formation efficiency (yU ^ 20 Gyr"') and an IMF 
flatter than in the solar neighbourhood. The MD they pre- 
dict peaks around solar metallicity and is shifted toward lower 
metallicities with a more extended metal poor tail with respect 
to our bulge observed MD. This disagreement was expected 
given the difference s between our bulg e MD and the one of 
Zoccah et all (I2q03h andlF ulbright et alJ (I2OO6I) fliat were used 
as references bv IBallero et al.l (l2007bh . In fact, to predict near 
solar metallicities MD consistent with high a-elements ratios, 
Ballero et alj (l2007bl) had to assu me an IMF flatter tha n in the 
solar neighbourhood b ut also ( see Balle ro et alj (|2007a)) flatter 
than the steUar IMF oflKroupa (200 1) and die IMF integrated 
over galaxies oflWeidner & Kroupa| (l2005l) . According to the 
tests done by the authors, to shift the predicted MD to more 
metal-rich (as observed in the present paper), the IMF should 
be further flattened. However, for IMF slope lower than 0.35, 



V. Hill et al.: The metallicity distribution of bulge clump giants in Baade's Window 



19 



cd 



^ 

Z 



25 
20 
15 
10 
5 




-,R^ 


^ 

>> 




N 

X 




\ 
\ 

N 
\ 


1 





0.5 






iac 



^ 



-0.5 



-0.4 -0.2 



0.2 

[Mg/Fe] 



0.4 0.6 




0.5 

[Mg/H] 



Fig. 15. Left panel: [Mg/Fe] distribution for the stars with very high probabilities (> 85%) to belong to population A (empty 
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the position of the peak of the predicted MD becomes insensi- 
tive to further slope decreases. Moreover, a flatter IMF would 
in turn lead to predict higher [Mg/Fe] and [O/Fe] ratios than 
observed. Thus, it seems there does not exist any appropriate 
set of fundamental pa rameters for the bulge chemical model of 
Ballero et al.l (l2007bl) to reproduce both our red clump MD and 
the [Mg/Fe] (or [O/Fe]) trend observed in the bulge. 

On the other hand, if we consider the observed MD as the 
result of two populations (see Se ct. I4.2ll, it would be m ore ap- 
propriate to compare the MD of iBallero et al.l (2007b) (given 
their bulge chemical model formation hypothesis) with the one 
associated to the old bulge component. This comparison is il- 
lustrated by the Fig. [16] where the old bulge population (A) 
is modelled by a Gaussian with a mean [Fe/H] of -0.3 and a 
width of 0.2 5 . With its actual parametrisation, the model of 
Ballero et al.l (2007b) also fails to reproduce the metal-poor 
population of the bulge. The predicted MD peaks at [Fe/H] val- 
ues higher than the mean of population A, and is wider than the 
Gaussian. It underpredicts the fraction of stars in population A 
at intermediate metallicities (-0.8 < [Fe/H] < -0.1) and pre- 
dicts a longer low metallicity tail ([Fe/H < -0. 8). However, 



based o n the parameter variations presented by B allero et al 
(l2007bl) . we would expect that it should be possible to find a 



set of parameters that would reproduce both the high [Mg/Fe] 
and metallicity distribution of population A. In particular, there 
would not be the need for a flatter IMF in the bulge anymore. 

5.2. Comparing to a chemodynamical model 

Figure[T6llower panel shows the comparison between our bulge 
iron MD a nd the predic t ed MD from two models (B and F 
models) of llmmeli et al.l (1200 41) convolved with the observa- 
tional errors. Ilmmeli et al.l (12004) investigated different scenar- 



phases (gas and stars) interstellar medium. The Galactic bulge 
is formed in very different way depending on values adopted 
for the efficiency of energy dissipation of the cold cloud com- 
ponent from which stars form. With a high efficiency of energy 
dissipation, a massive bulge is formed from the central merging 
of clumps of stars and gas at relatively early times (model B). 
On the contrary, with low efficiencies values the instability oc- 
curs in the stellar disc at later time forming a stellar bar which 
evolves and forms a pseudo-bulge (model F). 

The two predicted MD are completely different from the 
observed one. As already explained in the previous section, 
this was expected for the model B leading to the formation of a 
massive bulge. However, the B model MD does not agree either 
with the metal poor component of the observed MD (popula- 
tion A). It peaks around [Fe/H] = -0.40, close to the mean of 
the Gaussian representative of the old bulge but is wider than 
this population, and overproduces the fraction of metal-poor 
stars with respect to the actual observed MD. 

The observed and model F MDs peak at distinct metallic- 
ities ([Fe/H] — 1-0.3 and ~ -0.2 for the observed and model 
F respectively), but have a quite similar global shape. Indeed, 
if we artificially shift the model F MD of +0.5, the diff'erences 
between the two MD are less pronounced, even if differences in 
the lowest and highest metallicities regimes persist. The shifted 
F MD would overestimate the number of stars at low metallic- 
ities ([Fe/H] < -0.8 dex) and underestimate the number of 
stars at high metallicities ([Fe/H] > 0.1), in addition to show- 
ing a more progressive decrease of the number of stars than 
observed at high metallicities. We will come back in Sect l5.3l 
on these differences. However, if we believe the scenario of 



Immeli et al.l (120041) for the bulge formation (model F), we can 



ios of bulge, disc and halo populations formation from the evo- 
lution of a star-forming disc with a chemodynamical 3D model 
including initial conditions from ACDM results and a two 



not explain why such a large shift (+0.5) in mean metallicity 
would be required. Uncertainties in the stellar yields adopted 
can not explain such a large difference. To obtain such metal- 
rich stars, we would suggest that the gas forming the stellar disc 
(from which the pseudo-bulge is formed) could have been pre- 
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Fig. 16. Upper panel. Comparison between the bulge MP pre- 
dicted by the chemical model o f Ballero et al. (2007b) (red 



line) and by the model of llmmeli et al. (2004) (green line for 
the model B and blue hne for the model F ) with the observed 
MD (shaded histogram). The three model MDs have been con- 
volved with observational errors. Lower panel: . The previous 
model MDs were normalised to 50% and compared to the MD 
of the two subpopulations A and B (each of them correspond- 
ing to 50% of the total) illustrated here by two Gaussians with 
position and dispersion according to Tab. |5] 



viously enriched. Several factors could have contributed to the 
enrichment of this gas, but one possibility could be enriched 
gas left over (or blown out) from the old bulge formation. 

5.3. An attempt at modelling the two populations 

Without any satisfactory bulge formation model that could ex- 
plain the observed MD, we tried to build a simple modelling 
of the two populations according to t heir respective formatio n 
mechanism suggested by the results of lBabusiaux et al.l (120101) . 
The metal-poor population has kinematics that show no influ- 
ence of the Galactic bar, and in addition, we find it to be en- 
riched in [Mg/Fe], which reinforces out interpretation of this 
population as an old spheroid with a rapid time-scale forma- 
tion. To produce the chemical evolution of this population, we 



used a simple closed-box model. This model predicts that the 
metallicity (Z) distribution function (MDF) follows the relation 
f(Z) = p 'e"^''P where p is the stellar yield given by p =< Z >. 
For p, we adopted the mean values of the Gaussian associ- 
ated to the metal old population: -0.06 and -0.30 [Mg/H] and 
[Fe/H] respectively. The resulting [Fe/H] and [Mg/H] metal- 
licities distribution (scaled to 50% of the total population to 
represent population A only)) are compared to the red clump 
ones on Fig.[T7l( mid dle panels). 



As explained in iBabusiaux et al.l (120101) . the kinematical 



characteristics of the richest population suggest that it is un- 
der the influence of a bar. Adding to this the fact that [Mg/Fe] 
in this population is near solar, we suggest a formation over a 
long time scale through the evolution of the bar, itself originat- 
ing in disc instability. Within this scenario, population B stars 
observed in Baade's window would have been ejected from the 
inner regions of the Galactic disc. We then expect the chemi- 
cal composition of these stars to be the same as the one of the 
old inner disc at the time when the bar formed. However, as 
already mentioned, little is known about the inner regions of 
the Galactic disc. A metallicity gradient has been observed in 
young populations (OB stars, HII regions, Cepheids) but ob- 
servational constrains of the old populations are missing. To 
reproduce the MD of the richest population, we adopted the 
MD of the solar neighbourhood from Fuhrmann (2008), scaled 
it to 50 % of the total sample (to represent population B only), 
and shifted toward higher metallicity in o rder to r eproduce the 
possible gradient in the Galactic disc. The lFuhrm ann ( 2008) Fe 
and Mg metallicity distributions are symmetric around mean 
values of [Fe/H] = -0.03 + 0.01 and [Mg/H] = 0.00 ± 0.01, 
and narrower than the bulge MD (see Fi g.fTTl left panels). On 
the middle panels of Fig. [17] we plot the iFuhrmannI (l2008h Fe 
and Mg MD scaled to 50 % of the total bulge population, and 
shifted by +0.35 for both Fe and Mg, so that their mean val- 
ues coincide with those of the corresponding Gaussian com- 
ponents. These Mg and Fe MD shifted disc components seem 
to be good candidates to reproduce the observed bulge metal 
rich component. The main difference (more pronounced for Fe 
than for Mg MD) appears at very high metallicities, where the 
observed bulge MD shows a sharper cutoff than in that of the 
shifted disc. Assuming that this mock disc would originate in 
the very inner parts of the Galactic disc, around 8 k pc from the 



Sun, t he shifts applied to the local thin disc MD of lFuhrmann 
(2008) correspond to gradients of -0.047dex.kpc"' for Fe and 
-0.043dex.kpc"' for Mg. These gradients are similar to the gra- 
dients that have been obs erved in the youngest p opulations of 
the inner and outer disc ( IDaflon & Cunhall2004l) . However, in 
our case we are not interested in the value of the present day 
gradient in the Galactic disc but rather what it was in the past, 
when the stars that we now see in the pseudo-bulge would have 
formed. Recent observations of planetary nebulae a nd open 
clusters in the Galatic disc (see lMaciel & Costall2009l and ref- 
erences therein) suggest a flattening of the metallicity gradients 
with time which is not consistent with the gradient we derived. 
The chemical evolution models of the disc assuming an inside 



out f ormation scenario (iCescutti et alJ l2007t IChiappini et al 



200 lb predict gradients at the present day whose values agree 



with the observations of young and intermediate-age popula- 
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Fig. 17. The left panels show the comparison of the [Fe/H] M D (upper) and [Mg /H] MD (lower) with the thin (blue) and thick 
(green) discs MD in the solar neighbourhood from the study of iFuhrmannI (120081) . The distinction between the three populations 
is clear from their [Fe/H] and [Mg/H] MD in their global shapes, their mean values and the metallicity range they cover. In the 
middle panels the bulge [Mg/H] and [Fe/H] MD are compared to a close box model with a yield p - -0.27 and -0.06 for Fe and 
Mg respectively (black line) and to the previous thin disc MD shifted by +0.35 (blue) according to the results of the Gaussian 
separation (see Sect. 14.21 and Table |5]l. The combination of these two MD convolved are compared to the observed ones (right 
panels) for Fe (upper) and Mg (lower). For both elements, the model MD shows a very good agreement with the observed one. 
A small difference (compared to the relative errors) exists on the high metallicities tails: the high metallicity tail of the Mg and 
Fe MD decrease more sharply than the model distributions. 



tions of the disc. But, these mo dels also predict that the gra- 
dients were flatter at in the past. ICescutti et al.l (120071 private 
communication) found that 9 Gyrs ago, gradients were only of 
-0.014dex.kpc"' for Fe and of -0.017dex.kpc"' (in the range 
of galactocentric distances from 4 to 14 kpc). These values, if 
extended would also be inconsistent with our findings. 

On Fig. [17] we plotted the sum of the two modelled MD 
associated to the metal poor and metal rich components. The 
MD predicted by the Simple model have been convolved with 
the observed mean uncertaintie s for o ur sample (Fig. |6ll. The 
internal errors of the iFuhrmannI ("2008) data being much lower 
than the observed one, the disc MD have also been convolved 
with a Gaussian with cr = 0.18 forFe and cr = 0.12 for Mg. For 
both elements, the combined MD show a fairly good agreement 
with the observed ones. The main differences concern the metal 
poor tail and the richest part of the distributions. 

In the very metal poor regime ([Fe/H] < -0.7 and 
[Mg/Fe] < -0.5), the closed box model predicts more stars 
with respect to what is observed. The same lack is observed 



in the solar neighbourhood thin disc and is well known as the 
G dwarfs problem. For our red clump sample, it is not likely 
that this difference could come from a selection effect, at least 
down to [Fe/H] = -1.5 (see Sect. 12. 3t . As for the local disc, the 
G dwarf problem in the bulge could be partly or totally solved 
by including an infall of primordial gas in chemical evolution 
models. Additional observations of the bulge in the very metal 
poor regime would also help to answer this question. 

At higher metallicities, the observed MDs show a peak 
more pronounced and a sharper decrease than the shifted MD 
of the thin disc. These differences are quasi insensitive to 
changes up to 0.05 dex on the adopted shifts for Mg and Fe 
but could be reduced by decreasing the adopted width of the 
convolution (reproducing observed errors) down to an unreal- 
istic level of 0.10. In Sect. 14.1.31 we have argued that the high- 
metallicity tail of the iron MD may be biased at the highest 
metallicity stars, in a direction that is not yet very clear: ex- 
periments with synthetic spectra suggest that we may overes- 
timate the metallicities of super-solar metallicity stars by up 
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to 0.15 dex, whereas the estimated errors and sharpness of the 
decrease of the metal-rich tail may be read as a hint that the 
highest metallicities are underestimated. 

6. Summary and conclusions 

We have analysed a sample of 219 bulge red clump stars in 
Baade's window from R=20000 resolution spectra obtained 
with FLAMES/GIRAFFE at the VLT. This sample has been 
selected based on the location of stars in the Colour-Magnitude 
Diagram, in order to minimise the contamination by other 
Galactic components (< 10 % from the l atest p opulation 
synth esis Besangon model (Robi n et alJl2003uPicaud & Robin 
20041) . and hence is a good tracer of the bulge metallicity dis- 
tribution (MD) in the Baade's window. 

For these stars, with an automatic procedure, we have de- 
rived the stellar parameters and Fe abundances using an iron 
linelist established differentially to yuLeo and for a subsample 
of 162 stars (built by excluding stars with high uncertainties 
on [Fe/H]), we also derived the Mg abundances with spectral 
synthesis around the Mg I triplet at /I 63 19 A. 

The Fe and Mg metallicity distributions are both asym- 
metric, with mean values of +0.05 + 0.03 and +0.14 + 0.03, 
and median values of 0.16 and 0.21 respectively. They show 
a small proportion of stars at low metallicities (< 3 % with 
[Fe/H] < -0.7 or [Mg/H] < -0.4), only extending down to 
[Fe/H] = -1.1 or [Mg/H] = -0.7. The iron MD is not well 
reproduced by the predictions of the chemical mod els of the 
bulge available in the literature (IBallero et al.ll2007blla[) . neither 
by models adding dynamical aspects t o chemical prescriptions 
such as the one of jlmmeli et al.ll2004l) . 

The decomposition of the observed Fe and Mg MDs into 
Gaussian components has revealed the presence of two popula- 
tions of equal sizes (50% each): a metal-poor component cen- 
tred around [Fe/H] = -0.30 and [Mg/H] = -0.06 with a large 
dispersion and a narrow metal-rich component centred around 
[Fe/H] = 0.32 and [Mg/H] = +0.35. The metal poor compo- 
nent shows high [Mg/Fe] ratios (around 0.3) whereas stars in 
the metal rich c omponent are found to have near solar ratios. 
Babusiaux et al.l (|2010) also found kinematical differences be- 
tween the two components: the metal poor component shows 
kinematics compatible with an old spheroid whereas the metal 
rich component is consistent with a population supporting a 
bar In view of their chemical and kinematical properties, we 
suggest different formation scenarios for the two populations: 
a rapid formation timescale as an old spheroid for the metal 
poor component (old bulge) and for the metal rich component, 
a formation over a longer time scale driven by the evolution 
of the bar (pseudo-bulge). This latter component represents 
50% of the total Baad e's Window population, and is of old age 
(IClarkson et alJl2008l as indicated by deep colour-magnitude 
diagrams, see e.g.), indicating that a fairly massive disc was 
present some lOGyr ago, which may turn out to be a strong 
constraint on models of disc formation. It is suggestive that the 
local thick disc and old bulge display the same [Mg/Fe] en- 
hancement, pointing at similarly rapid star formation process. 
At this point, it is however not possible to constrain further the 
relationship that the old bulge bears with the thick disc, both 



because of lack of constraints on the central parts of the thick 
disc (both chemistry and kinematics), and because our bulge 
sample being located on the minor axis, we have a very poor 
constraint on rotation and the disky nature of this component 
can therefore not be probed from our data. We however not that 
in their recent study, iBensbv et alJ (l2010al) finds that the inner 
parts of the thick disc is chemically very similar to that in the 
solar neighbourhood, lending strong support to the idea that the 
metal-poor bulge is indeed chemically related to the thick disc. 
We also note that the numerous metal-rich globular clusters of 
the Milky Way, that are concentrated towards the inner parts of 
the galaxy, would be classified as old bulge in our scheme (see 
e.g. Bica et al. 2006 that show clusters with [Fe/H] < -0.75 
are confined within 2kpc of the centre). It is interesting to note 
that the most metal-rich globular cluster known so far is NGC 
6528 (Ortolani et al. 2007), with [Fe/H] = -0.2 (ZoccaU et al. 
2004), therefore within the metal-poor bulge component iden- 
tified here. 

Guided by these results, we built a simple model combin- 
ing two components. We used a simple closed box model to 
predict the metal poor population contribution, whereas we rep- 
rese nted the metal ric h population by the MD of the local thin 
disc JFuhrmannlbOOSh . shifted to the centre of the correspond- 
ing metal rich population. The resulting combined MD shows a 
good agreement with the red clump one on the whole metallic- 
ity range except for the very metal poor and very metal rich 
regimes, where more stars are predicted than observed. The 
bulge (old bulge) thus seems to exhibit the G dwarfs problem. 
The shift we had to apply to the local thin disc MD to fit the 
metal-rich population is however quite large (although not im- 
possible) when compared with plausible values of the gradients 
in the early Galactic disc. One possibility could be that these 
stars were born in a pre-enriched environment, enriched by the 
remnant of the old bulge. 

The scenario(s) proposed here should now be turned into 
proper modelling (chemical and kinematical) of the bulge that 
would include two formations events (linked or not), to eval- 
uate their viability. Expanding similar combined chemical and 
kinematical studies of large samples towards various directions 
in the Galactic bulge would also be of great help for answering 
these questions. In particular, probing the bulge at lower lati- 
tudes and along longitude should trace the relative importance 
of the metal-rich and metal-poor populations and definitively 
establish their proposed different (kinematical and chemical) 
nature as an old bulge and pseudo-bulge. Because these re- 
gions are much more heavily obscured than Baade's window, 
this kind of observations will have to be carried out in the in- 
frared. Finally, we note that our group has undertaken similar 
red clump observations of the inner parts of the Galactic disc 
that are expected to shed new light on the true relation between 
the disc and the pseudo-bulge. 
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